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1 Introduction 

In this paper we study concentration phenomena for the problem 

-e 2 Au + u = u p in fL 



(Pe) 



|^ = on dQ, 

u > in $7, 



where Q. is a smooth bounded domain of R , p > 1, and where v denotes the unit normal to 
<9$7. Given a smooth embedded non-degenerate minimal submanifold K of d£l, of dimension 



A; G {1, ... ,N — 2}, we prove existence of solutions of (i^J concentrating along i^. Since the 
solutions we find have a specific asymptotic profile, which is described below, a natural restriction 
on p is imposed, depending on the dimension N and k, namely p < A ~ ; +~ 



N-k-2- 



Problem ( |P e [ ) or some of its variants (including the presence of non-homogeneous terms, 
different boundary conditions, etc.) arise in several contexts, as the Nonlinear Schrodinger 
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Equation or from modeling reaction-diffusion systems, see for example [3], [22], [46] and refer- 
ences therein. A typical phenomenon one observes is the existence of solutions which are sharply 
concentrated near some subsets of their domain. 

Concerning reaction-diffusion systems, this phenomenon is related to the so-called Turing's 
instability, [55] . According to this principle, reaction-diffusion systems whose reactants have 
very different diffusivities might generate stable non-trivial patterns. This is indeed more likely 
to happen when more reactants are present since, as shown in [12], [42], scalar reaction-diffusion 
equations in a convex domain admit only constant stable equilibria. 

A well-know system is the following one 



(GM) 



U t = d 1 AU-U+^ inftx(0,+oo), 
V t = d 2 AV-V + ^ inOx(0,+oo), 
M = ™=0 ondQx (0,+oo) 



introduced in [25] to describe some biological experiment. The functions IA and V represent 
the densities of some chemical substances, the numbers p, q, r, s are non-negative and such that 
< < an d it is assumed that the diffusivities d\ and d 2 satisfy d\ <C 1 <C d 2 . In 
the stationary case of \GM\ . as explained in [46], [49], when d 2 — > +oo the function V is close 
to a constant (being nearly harmonic and with zero normal derivative at the boundary), and 



therefore the equation satisfied by U is similar to ( P £ ) , with e = d\ 



The typical concentration behavior of solutions u £ to QP £ D is via a scaling of the variables 
in the form u e (x) ~ no (^f^J > where Q is some point of ft, and where uq is a solution of the 
problem 

(1) -Au + u = u p Q inR N (or in = {(xi,..., x N ) € M N : x N > 0}), 

the domain depending on whether Q lies in the interior of f2 or at the boundary; in the latter 
case Neumann conditions are imposed. 

When p < -j^l (and indeed only if this inequality is satisfied), problem ([I]) admits positive 
radial solutions which decay to zero at infinity. Solutions of ( |P £ [ ) with this profile are called spike- 
layers, since they are highly concentrated near some point of £1. There is an extensive literature 
regarding this type of solutions, beginning from the papers [35], [37], [38]. Indeed their structure 
is very rich, and there are also solutions with multiple peaks, both at the boundary and at the 
interior of O. We refer for example to the papers [II], [U], [26], [27], [28], [29], [33], [3J], [58]. 



In recent years, some new types of solutions have been constructed: they indeed concentrate 
at sets of positive dimension and their profile consists of solutions of (pQ) which do not decay to 
zero at infinity. In [39], [10] it has been shown that given any smooth bounded domain 17 C M. N , 
N > 2, and any p > 1, there exists a sequence Sj — > such that {P Ej ) possesses solutions 
concentrating at dfl along this sequence. Their profile is a solution of (pQ) (for = 1) on the 
half real line which tends to zero at infinity and which satisfies the condition u' (0) = 0. This 
function can also be trivially extended as a cylindrical solution to ([T]) on the whole . 

Later in [38] it has been proved that, if is a smooth bounded set of M 3 , if p > 1 and if 
h is a closed, simple non-degenerate geodesic on d£l, then there exists again a sequence (£j)j 
converging to zero such that (P £j ) admits solutions u £j concentrating along h as j tends to 
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infinity. In this case the profile of u £ . is a decaying solution of (pQ) in again extended to a 
cylindrical solution in higher dimension. 

These are examples of a phenomenon which has been conjectured to hold in more general 
cases: in fact it is expected that, under generic assumptions, if Q C R N and if k is an integer 



between 1 and N — 1, there exist solutions of (P £ ) concentrating along /c-dimensional sets when 



e tends to zero. While the case k = N — 1 has been tackled in [40] , the goal of the present paper 
is to consider k < N — 2, and to prove this conjecture under rather mild assumptions on the 
limit set. Before stating our main theorem we introduce some preliminary notation. 

Given a smooth A;-dimensional manifold K of dQ,, and given any q £ K we can choose a 
system of coordinates (y, Q) in £1 orthonormal at q and such that (y, 0) are coordinates on K, 
and with the property that 

d d d 
(2) -r—\q&T q K, o = l,...,fe; — \ q eT q dtt, i = l,...,n; — \ q = v(q), 

where we have set n = N — k — 1 . Our main theorem is the following: we refer to Section [2] for 
the geometric terminology. 

Theorem 1.1 Let SI C Mr, N > 3, be a smooth and bounded domain, and let K C Oft be a 
compact embedded non- degenerate minimal submanifold of dimension k £ {1, ... , N — 2}. Then, 

if p G (l, jvZjfej^ ) ' there exists a sequence Ej — > such that (P £j ) admits positive solutions u e . 
concentrating along K as j — ► oo. Precisely there exists a positive constant C, depending on 

dist{x ,K) 

f2, K and p such that for any u £j {x) < Ce Ce J ; moreover for any q 6 in a system 

>rdin> 

wo : — > M is the unique radial solution of 



of coordinates (y,() satisfying ([2]), for any integer m one has u Ej (0,£j-) loc — t wo{'), where 



' -Au + u = U p inR n+1 

(3) 



^ = ondRl +1 , 



u > 0,u £ ^ 1 (R" +1 ). 



Remarks 1.2 (a) Differently from the previous papers concerning the case N = 3 and k = 1, 
or concentration at the whole d£l, we require an upper bound on p depending on N and k. This 
condition is rather natural, since © is solvable if and only if p < fEfjf , see fM/, J^, Jjgjj 
and in this case the solution is radial and unique (up to a translation), see ]23$ , J22F- In any 
case, our assumptions allow supercritical exponents as well. 

(b) As for the results in JgffJ /. [39] and 1401, existence is proved only along a sequence Ej — * 
(actually with our proof it can be obtained for e in a sequence of intervals (a,j,bj) approaching 
zero, but not for any small e). This is caused by a resonance phenomenon we are going to 
discuss below, explaining the ideas of the proof. This resonance is peculiar of multidimensional 
spike-layers, see also [20], and other geometric problems, see \31^ , \43^ . In some cases, when 
some symmetry is present, it is possible to get rid of this resonance phenomenon working in 
spaces of invariant functions. We refer for example to the papers Q 0, 0, [13 [IB [73 \4~°1 - 
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We can describe the resonance phenomenon, which causes the main difficulty in proving Theorem 
1 . H in the following way. By the change of variables x *— > ex, we are reduced to consider the 
problem 



—Au + u = vP in Q E , 
|^=0 ondn 
u > in 



ei 



where VL e = -O. As for ((2J), given q G K £ := ~K, we can choose scaled coordinates (y, C) 
on S\ such that <9 y Jq G T^K £ , d^\q G T^d£l £ and e?f n+1 |$ = ^(g). Then, letting tt £ denote 
the scaling of u e to fi £ , we have that, in a plane through q normal to K E , u £ behaves like 
u £ (Q,() = u £ (0,e() ~ ifo(C)- This amounts to the fact that ^ u>o (dist(x, K £ )), x G fi e , 

and therefore u e has a fixed profile in the directions perpendicular to the expanding do ma in K £ . 



Since the function wq (dist(x, K £ )) can be considered as an approximate solution to (P E ), it is 



natural to use local inversion arguments near this function in order to find tr ue so lutions. For this 



purpose it is necessary to understand the spectrum of the linearization of ( P e 1 at approximate 
solutions. 

For simplicity, let us assume for the moment that K is (N — 2) -dimensional, namely that its 
codimension in d£l is equal to 1, as in [38]. Then, letting v denote the normal to K in d£l, we 
can parameterize naturally a neighborhood of K £ as a product of the form K £ x (— | , |), where 
5 is a small positive number, via the exponential map in d£l £ 



(4) (y, s) ^ exp^(sz>); (y, s) G K £ x 



6 6_ 



Similarly, if v (y, s) is the inner unit normal to d£l £ at the image of (y, s) under the above map, 
we can parameterize a neighborhood of K £ in 0, £ with a product K £ x (— |, |) x (0, |) by 

(y,s,t) ^exp° QB (sv)+tv(y,s); (y,s,t)£K £ x x f°> Q • 

When e tends to zero, the standard Euclidean metric of f2 £ becomes closer and closer (on the 
above set) to the product of the metric of K £ and that of M 2 (parameterized by the variables 
s and t as cartesian coordinates). Therefore, since the set (— f , f) x (0, |) converges to = 



{(s,t) G R 2 : t > 0}, in a first approximation we get that the linearization of (|P g [) at u £ is 



(5) 



-Ak e u — d 2 s u — d 2 t u + u — pwo(()u = in K £ x 
|^ = on K £ x 



The spectrum of this linear operator can be evaluated almost explicitly. Referring to Section U] 
for details (see also [38], Proposition 2.9 for the case N = 3), here we just give some qualitative 
description of its properties. 

Given an arbitrary function u G H 1 (K £ x R?_), we can decompose it in Fourier modes in the 
variables £ = (s,t) as 

j 

Here 0j are the eigenfunctions of the Laplace-Beltrami operator on K, namely —Ax4>j = Pj4>ji 
j = 0, 1, 2, ... , where the eigenvalues (pj)j are counted with their multiplicities. 
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If u is an eigenfunction (with respect to the duality induced by the space H 1 ^^. x of 
the linear operator in ([5]) with corresponding eigenvalue A, then it can be shown (see Section H] 
for details) that the functions Uj satisfy the equation 

, R v f (1 - A) [-A Uj + (1 + a)uj] - pw^Uj = in M.%, 

[ j 1^ = on dR 2 + , 

where a = £ 2 pj- It is known that when a = the latter problem admits a negative eigenvalue 
770 (with eigenfunction wq), a zero eigenvalue <jq (with eigenfunction d s wo), while all the other 
eigenvalues are positive. This structure is due to the fact that u>o is a mountain-pass solution of 
([3l) (so its Morse index is at most 1), and the presence of a kernel derives from the fact that this 
equation is invariant by translation in the s variable. When a is positive instead, it turns out 
that the first eigenvalue rj a of ([6]) and the second one a a are strictly increasing functions of a 
with positive derivative, and tend to 1 as a — > +00; moreover, the eigenfunctions corresponding 
to r] a (resp. a a ) are radial (resp. odd in s) for every value of a. In particular, there exists a > 
such that rfa = 0, so when e 2 pj is close to a we obtain some small eigenvalues of the original 
linearized problem (|5"|) . 

From the monotonicity in a and from the Weyl's asymptotic formula for pj, it follows that 

2 

the eigenvalues of the operator in (J3J) are, roughly, either of the form r]o + e 2 j N ~ 2 for some j G N, 

9 2 

or of the form e l N ~ 2 for some I G N, or have a uniform positive bound from below. 

In the case of general codimension it is not possible to decompose a neighborhood of K 
(in <9fi) as for (j4]), but instead one has to model it on the normal bundle of K £ in U £ , see 
Subsection 14.21 for details. Considering the corresponding approximate linearized operator, one 
can prove that its eigenvalues are now, roughly either of the form i] £ 2 p . ~ 770 + £ 2 j k , or of the 

form & e 2 UJl ~ e 2 /fc, j, I G N, or, again, have a uniform positive bound from below. Here (pj)j 
still represent the eigenvalues of the Laplace-Beltrami operator on K, while the numbers {loi)i 
stand for the eigenvalues of the normal Laplacian of K (considered as a submanifold of dfl), see 
Section [2] for its definition and the corresponding Weyl's asymptotic formula. We are interested 
in particular in the following two features of the spectrum: 

1) resonances: there are two kinds of eigenvalues which can approach zero. First of 

„ 2 

all, those of the form n a when a is close to a. This happens when e j k — a, namely when 

j ~ e~ k ; furthermore, the average distance between two consecutive such eigenvalues is of order 

e 2 jk~ l ~ j~ l ~ e k . The other resonant eigenvalues are of the form a a ~ a for a close to zero, 

namely when a = e 2 lk and I is sufficiently small (compared to, say, some negative power of e). 

Hence the distance from zero of the smallest eigenvalues of this type is of order e 2 . Indeed, an 

accurate expansion in e, see Subsection 15.21 yields that this distance is bounded from below by 

a multiple of e 2 when K is a non-degenerate minimal submanifold. 

2) eigenfunctions: as for the case of codimension 1, it turns out that the eigenfunctions 
corresponding to the 77 a 's are of the form (pj(ey)uj(Q, where Uj is radial in the variable £ (£ 
represent here some orthonormal coordinates in the normal bundle of K £ ). The function (pj 
instead oscillates faster and faster as e tends to zero, since j is of order e~ k . On the other hand 
it is possible to show, see Subsection 14.21 that the eigenfunctions corresponding to the a a 's are 
products vi(\C\)(Ci^Pi)n, where {-,-)n is the scalar product in NK, and where (pi is a section of 
the normal bundle NK e , and precisely an eigenfunction (scaled in e) of the normal Laplacian 
of K. Since the resonant modes correspond to low indices I, (pi does not oscillate as fast as the 
resonant 0j's. 
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So far we considered an approximate operator, because in ([5]) we assumed a splitting of the 
metric into a product. Since we expect to deal with small eigenvalues, a careful analysis of 
the approximate solutions is needed (to apply local inversion arguments), and also a refined 
understanding of the small eigenvalues with the corresponding eigenfunctions. 

Therefore we first try to obtain approximate solutions as accurate as possible. For doing 
this , as in \38 \ \39 \ 140] . one can introduce suitable coordinates on £l £ near K £ , expand formally 



\P £ \ in powers of e, and solve it term by term using functions of the form 
(7) 

ui,e(y,()= [w + ew 1 + --- + e I w I ]{ey,C' + ^o{ey) + --- + e I 2 $/- 2 (ey), Cn+i); C=(C',Cn+l)- 

Here $0, • • • , <E>/_2 represent smooth sections of the normal bundle NK, and the functions (wi)i 
are determined implicitly via equations of the type 



(8) 



-Awi + Wi -pw (()wi = Fi(ey,w , ■ ■ ■ ,11^-1, $ , • • • ,^-2) in R" +1 , 
^ = on dRl +1 . 



Notice that the operator acting on Wi is nothing but the linearization of ([3]) at wq (shifted 
in by <3?o + " " " + e 2 ^7— 2); which has an n-dimensional kernel due to the invariance by 
translation in The functions 3>i are chosen in order to obtain orthogonality of F{ to the 
kernel, and to guarantee solvability in Wi. In doing this, the non-degeneracy condition on K 
comes into play, since the <3>i's solve equations of the form Z^i = Gi(y). Z denotes the Jacobi 
operator of K, related to the second variation of the volume functional, which is invertible by 
the non-degeneracy assumption on the minimal submanifold. Notice also that we wrote the 
variable y with a factor e on the front. This is in order to emphasize the slow dependence in y of 
these functions. In fact, recalling that (in the model problem described above) resonance occurs 
mostly when dealing with highly oscill atin g eigenfunctions, if we require slow dependence in y 



then there is no obstru ction in solving (^) up to an arbitrary order e 1 . 

Next one linearizes (Pr_ ) near the approximate solutions just found. Compared to the above 
model problem, the eigenvalues will be perturbed by some amount, due to the presence of the 
corrections (u)i)i and to the geometry of the problem. In fact the amount will be in general 
of order e, since this is the size of the corrections (from the w^s and the expansions of the 
metric coefficients, see Lemma \3.2\i . This prevents a direct control of the small eigenvalues of 
the linearized operator (at uj t£ ) since, as discussed above, the characteristic size of the spectral 
gaps at resonance are of order e 2 or e k . 

To overcome this problem, we look at the eigenvalues as functions of e. The counterparts of 
the numbers o~ £ 2 Wl can be again obtained via a Taylor's expansion in e, and they turn out to be 
constant multiples of e 2 times the eigenvalues of Z (up to an error of order o(e 2 )), so they are 
never zero. On the other hand, the counterparts of the ^ e 2 pj 's could vanish for some values of e 

9 - 

but, recalling the expansion rj £ 2 p . ~ r/Q + e j * , one can hope that generically in e none of these 
eigenvalues will be zero. 

This is indeed shown using a classical theorem due to T. Kato, see [30], pag. 445, which allows 
us to estimate the derivatives of the eigenvalues with respect to e. To apply this result one needs 
some control not only on the initial eigenvalues but also on the corresponding eigenfunctions, 
and this is what basically the last sections are devoted to. There we prove that if A = o(e 2 ) 
is an eigenvalue of the linearized operator, the eigenfunctions (up to a small error) are linear 
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combinations of products like (j)j(ey)uj((), for j ~ e~ k and for suitable functions Uj radial in £. 
Then we deduce that ^ is close to a number depending on e, N,p and K only. As a consequence, 
the spectral gaps near zero will shift, as e varies, almost without squeezing, yielding invertibility 
for suitable values of the parameter. This method also provides estimates on the norm of the 
inverse operator, which blows-up with rate max{e _fe , £~ 2 } when e tends to zero, see Remark 16.81 
Finally, a straightforward application of t he im plicit function theorem gives the desired result. 



To fix the ideas, when p < ]^r§, solutions of ( |P E [ ) can be found as critical points of the following 
functional 

(9) Je(u) = \( (\Vu\ 2 + u 2 )-^— I K+\ u€H\Q £ ). 

rife 

One proves that || J' £ (ui t£ )\\jfirQ e \ < Cj^e a for e small. Even if the norm of the inverse 
linear operator blows-up when e tends to zero, choosing I sufficiently large (depending only on 
k and p), one can find a solution using the contraction mapping theorem near ui e . 

The general strategy of this proof, and especially Kato's theorem, has been used in [38], 
[39j and |3Q], so throughout the paper we will be sketchy in the parts where simple adaptations 
apply. However the present setting requires some new ingredients: we are going to explain next 
what are the differences with respect to these and to some other related papers. First of all, 
compared to [39], [3D], where the case k = N — 1 was treated, here we need to characterize the 
limit set among all the possible ones, since the codimension is higher, and this reflects in the 
fact that the limit problem ([3]) is degenerate. This requires to introduce the normal sections 
<3?o, • • • j $1-2 m 0) and to use the non-degeneracy condition on K. 

The localization of the limit set has been indeed also faced in [38]. Here, apart from in- 
cluding that result as a particular case, allowing higher dimensions and codimensions, we need 
a more geometric approach. The main issue, as we already remarked, is that we cannot use 
parameterizations with product sets as in (jJJ, since the normal bundle of K is not trivial in 
general. At this point some interplay between the analytic and geometric features of the prob- 
lem is needed. In particular the first and second eigenfunctions of the linearization of ([3|) (the 
profile of u £ at every point q of K) can be viewed of scalar or vectorial nature. More precisely, 
the eigenfunction corresponding to the first eigenvalue is radial and unique up to a scalar mul- 
tiple. On the other hand the eigenfunctions corresponding to the second eigenvalue have the 
symmetry of the first spherical harmonics in the unit sphere of N q K, and they are in one-to-one 
correspondence with the vectors of N q K. The same holds true for the eigenfunctions of problem 
(J6j) when a > 0. When q varies over the limit set, these eigenfunctions (which are the resonant 
ones), depending on their symmetry determine respectively a scalar function on K or a section 
of the normal bundle NK, on which the Laplace-Beltrami operator or the normal Laplacian act 
naturally, see in particular Section 0J Apart from these considerations some other difficulties 
arise, more technical in nature, due to the more general character of the present result compared 
to that in [38]. Heavier computations are involved, especially since the curvature tensors have 
more components, and some extra terms appear. Anyway, some of the arguments have been 
simplified. 

Finally, we should point out the differences with respect to the papers [20] , [37], [43], where 
also special solutions of the Nonlinear Schrodinger equation or constant mean curvature surfaces 
are found. In ]20j and [43] the spectral gaps are relatively big, and the eigenvalues can be located 
using direct comparison arguments, so there is no need to invoke Kato's theorem. In [37] 
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arbitrarily small spectral gaps are allowed, but while there one has to study a partial differential 
equation on a surface only, here we need to analyze the equation on the whole space, which 
takes some extra work. Also, the Riemannian manifold we consider here, d£l, has an extrinsic 
curvature as a subset of R , and therefore some error terms turn out to be of order e, and not 
e 2 , see Remark 13.41 (a). Nevertheless, we take great advantage of the geometric construction in 
|37j . especially in their choice of coordinates near the limit set. We believe that our method 
could adapt to study concentration at general manifolds for the Nonlinear Schrodinger equation 
as well, as conjectured in [5J. 

The paper is organized in the following way. We first introduce some notations and conven- 
tions. In Section [2] we collect some notions in differential geometry, like the Fermi coordinates 
near a minimal submanifold, the normal Laplacian, the Laplace-Beltrami and the Jacobi oper- 
ators as well as the asymptotics of their eigenvalues. In Section [3] we construct the approximate 
solution ui t£ . In Section H] we study some spectral properties for the limit problem ([3]) (with 
some extension) and we then derive a model for the linearized operator at uj jE . In Section [5] 
we turn then to the real linearized operator: we construct some approximate eigenfunctions 
which allow us to split our functional space as direct sum of subspaces for which the linearized 
operator is almost diagonal. In Section [6J using this splitting we characterize the eigenfunctions 
corresponding to resonant eigenvalues. From these estimates we can obtain invertibility, via 
Kato's theorem, and prove our main result Theorem 11.11 
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Notation and conventions 

- Dealing with coordinates, Greek letters like a, /?,... , will denote indices varying between 1 
and N — 1, while capital letters like A, B, . . . will vary between 1 and N; Roman letters like a 
or b will run from 1 to k, while indices like . . . will run between 1 and n := N — k — 1. 

- {j, . . . , Q n , Cn+i w ih denote coordinates in M n+1 = M. N ~ k , and they will also be written as 

C' = (Cl,...,Cn), C=(C',Cn+l). 

- The manifold K will be parameterized with coordinates y = . . . , y k ). Its dilation K e := -K 
will be parameterized by coordinates (yi, . . . , yk) related to the y's simply by y = ey. 

- Derivatives with respect to the variables y, y or £ will be denoted by dy, d y , d^, and for brevity 
sometimes we might use the symbols da and <9j for dy a and 3q respectively. 

- In a local system of coordinates, (~g a g)ai3 are the components of the metric on dtl naturally 
induced by M. N . Similarly, (g^g )ab are the entries of the metric on in a neighborhood of the 
boundary. {H a p) a p will denote the components of the mean curvature operator of dQ into R . 
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Below, for simplicity, the constant C is allowed to vary from one formula to another, also within 
the same line, and will assume larger and lager values. It is always understood that C depends 
on Q, the dimension N and the exponent p. It will be explicitly written C[, C$, • • • , if the 
constant C depends also on other quantities, like an integer I, a parameter 5, etc. Similarly, the 
positive constant 7 will assume smaller and smaller values. 

For a real positive variable r and an integer m, 0(r m ) (resp. o(r m )) will denote a function 



for which 



also write o e (l) 



Q(r m ) 



remains bounded (resp. 



o(r m ) 



tends to zero) when r tends to zero. We might 



: : or a quantity which tends to zero as e tends to zero. With 0(r m ) we denote 
functions which depend on the above variables (y, C)> which are of order r m , and whose partial 
derivatives of any order, with respect to the vector fields d a , rdi, are bounded by a constant 
times r m . 

Li will stand in general for a differential operator of order at most i in both the variables y 
and C (unless differently specified), whose coefficients are assumed to be smooth in y. 

For summations, we might use the notation to indicate that the sum is taken over an 
integer index varying from [c] to [d] (the integer parts of c and d respectively). We might use 
the same convention when we make an integer index vary between c and d. We also use the 
standard convention of summing terms where repeated indices appear. 

We will assume throughout the paper that the exponent p is at most critical, namely that 
P < W^5> so that problem (i^) is variational in H 1 ^). We will indicate at the end what are 
the arguments necessary to deal with the general case. 



2 Geometric background 

In this section we list some preliminary notions in differential geometry. First of all we introduce 
Fermi coordinates near a submanifold of d£l, recall the definition of minimal submanifold, and 
introduce the Laplace-Beltrami and the Jacobi operators, together with some of their spectral 
properties. We refer for example to [6] and [53] as basic references in differential geometry. 

2.1 Fermi coordinates on d£l near K 

Let K be a /c-dimensional submanifold of (d£l, ~g) (1 < k < N — 1) and set n = N — k — 1 (see our 
notation). We choose along K a local orthonormal frame field ((E a )a=l,—k, {Ei)i=i, — ,n) which is 
oriented. At points of K, TOO, splits naturally as TK © NK, where TK is the tangent space to 
K and NK represents the normal bundle, which are spanned respectively by {E a ) a and {Ej)j. 

Denote by V the connection induced by the metric g and by V N the corresponding normal 
connection on the normal bundle. Given q S K, we use some geodesic coordinates y centered 
at q. We also assume that at q the normal vectors (-Ej)j, i = 1, ... ,n, are transported parallely 
(with respect to V^) through geodesies from q, so in particular 

(10) g(V Ea Ej ,Ei) =0 at q, i,j = l,...,n,a = l,...,k. 
In a neighborhood of q, we choose Fermi coordinates (y, Q on d£l defined by 

n 

(11) (y,C) -^expf (J2 &Ei); (y,C) = ((y a )a,((i)i), 

1=1 
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where exp^ is the exponential map at y in dU. 

By our choice of coordinates, on K the metric g splits in the following way 

(12) g(q) = g ab (q) dy a ®dy b + g^q) dQ <g> dQ; q G K. 
We denote by T^(-) the 1-forms defined on the normal bundle of K by 

(13) T b a (E i )=g(V Ea E b ,E l ). 

We will also denote by Rap^s the components of the curvature tensor with lowered indices, which 
are obtained by means of the usual ones Rp^$ by 

Ra/3j5 = 9 a <j R%f8- 

When we consider the metric coefficients in a neighborhood of K, we obtain a deviation from 
formula (|12p . which is expressed by the next lemma, see Proposition 2.1 in [37j for the proof. 
Denote by r the distance function from K. 

Lemma 2.1 In the above coordinates (y,C)> f or an U a = 1, and any i,j = 1, ...,n, we have 
^•(0,0 =5 ij + lR istj ( s t t + 0(r 3 ); 
g aj (0,C) =0(r*); 

9ab(0, C) = Sab ~ 2^) C, + [Rsabi + K(E S ) T^E^] ( S Q + 0(r 3 ). 
Here Ristj are computed at the point q of K parameterized by (0,0). 

2.2 Normal Laplacian, Laplace-Beltrami and Jacobi operators 

In this subsection we recall some basic definitions and spectral properties of differential operators 
associated to minimal submanifolds. We first recall some notions about the Laplace-Beltrami 
operator, the normal connection and the normal Laplacian. 

If (M, g) is an m-dimensional Riemannian manifold, the Laplace-Beltrami operator on M is 
defined in local coordinates by 

(14) A g = — L= d A ( ^/fetg~g AB d B ), 

where the indices A and B runs in 1, . . . , to, and where g AB denote the components of the inverse 
of the matrix gAB- 

Let K C M be a /c-dimensional submanifold, k < to — 1. The normal connection V w on a 
normal vector field V is defined as the projection of the connection onto NK. Moreover, 
one has the following formula regarding the horizontal derivative of the product (-,-)n in the 
normal bundle (see [53], Volume 4, Chapter 7.C, for further details) 

X(V, W) N = (V#V, W) + (V, V%W), 

for any smooth sections V and W in NK. If we choose an orthonormal frame (Ei)i for NK 
along K, we can write 

VlE j = (3Uc\)E l , 
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for some differential forms &\ (we recall our notation da = Sr)- Since the normal fields (Ei)i 
are chosen to be orthonormal, it follows that for any horizontal vector field X there holds 
X{Ei,Ej)N = 0, and hence one has 

(15) $ (8a) = (da) Vl,j = l,...,n:=m-k. 

This holds true, in particular, if we choose Fermi coordinates. Since indeed the normal fields are 
extended via (normal) parallel transport from q to some neighborhood through the exponential 
map, it follows that (3j(da)(0, 0, . . . , y a , 0, . . . , 0) =0, and hence 

(16) tf-(da)=0 &tq Va= l,...,k, andV/,j = l,...,n; 

(17) da(pUda)) =0 at q V a = 1, . . . , k, and V l,j = 1, . . . , n. 



Recalling these facts, we can derive the expression of the normal Laplacian in Fermi coordinates 
in the following way: given a normal vector field V = V 3 Ej, there holds 

VgV = d s ViE j + vW } (d s )E l . 

For any two normal vector fields V and W we have, by the definition of A^ 

f {v N vy N w) N dv-g = - f (A%v,w) N dv w . 

Jk Jk 
We compute now the expression of A^ evaluating the left-hand side and integrating by parts 

[ {V N V,V N W) N dV g - = [ (daV 3 E j + V 3 f3 l J (da)E h d E W i E i + W l ^(d E )) g ab V^V 
Jk Jk n ' n 

= I \daV i d I ;W i + daV 3 W i p{(d E )+V^ l j (da) ^W' 

Jk 1 

+ vm i {3) (da) ft (%)] g ab Vtet]j 
This quantity, for any V and W, has to coincide with — J K (A^-VyW l \/det g, so we deduce that 

(A K vy = A K (V') + -^d- b {y 3 f3)(d E )g ab V^ 



(18) - T b {daV j (3i (%) + (da) Pi (<%) ) v 7 ^- 

In Fermi coordinates at q, which is parameterized by (0,0), we have that 



(19) g ab = 5 a b, ckg a b = and V det 9 = 0, 

and we also have (|16|) - (|17| ). Hence the last formula simplifies in the following way 

(20) (A K VY = A K (V { ) at q. 
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Let C°°(NK) be the space of smooth normal vector fields on K. For $ G C°°(NK), we can 
define the one-parameter family of submanifolds t \— > Kt.§ by 

(21) K ti * := {expf (i$(y)) : y G K}. 

The first variation formula of the volume is the equation 



Vol(iT t)# )= / (^h^dVk-, 

t=o 



where h is the mean curvature (vector) of K in dCl, (•, -)jv denotes the restriction of g to iV-fT, 
and dVx the volume element of K. 

The submanifold i<C is said to be minimal if it is a critical point for the volume functional, 
namely if 



/ x d 

< 23 > * 



Vol(Kt,*) = for any $ G C°°\NK) 

t=0 



or, equivalently by (|22|) . if the mean curvature h is identically zero on K. It is possible to prove 
that, if T h a (Ei) is as in (JTHJ), then 

(24) K is minimal 44> Ta(i^) = for any i = l,...n. 

We point out that in the last formula we are summing over the index a, which is repeated. 

The Jacobi operator 3 appears in the expression of the second variation of the volume functional 
for a minimal submanifold K 



< 25 > £ 



Vol(K t ^) = - / (3$, <S>) N dV K ; $ G C°°(NK), 

t=0 JK 



and is given by 

(26) 3$ := -A$$ + 9^$ - 58^$, 

where SK^, 03^ : NK -> iVif are defined as 

= (R(E a , $)E a ) N ; g(<8 N $,n K ) := T a b ^)T b a (n K ), 

for any unit normal vector to K. The operator is the normal Laplacian on if defined 
in flUD. 

A submanifold K is said to be non- degenerate if the Jacobi operator 3 is invertible, or equivalently 
if the equation 3$ = has only the trivial solution among the sections in NK. 



We recall now some Weyl asymptotic formulas, referring for example to [13], or to [32] and 
for further details. Let (M,g) be a compact closed Riemannian manifold of dimension m, and 
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let A g be the Laplace-Beltrami operator. Letting (pi)i, i = 0, 1, ... , denote the eigenvalues of 
— A g (ordered to be non-decreasing in i and counted with their multiplicity), we have that 



(27) p ^ Cm {Vd(M)) 

where Vol(M) is the volume of (M,g) and C m is a constant depending only on the dimension 
m (the Weyl constant). A similar estimate, which can be proved using (|18p and (|27p . holds 
for the normal Laplacian A^ on a fc-dimensional submanifold K C M. In fact, letting (ujj)j, 
j = 0, 1, ... , denote the eigenvalues of — A^ (still chosen to be non-decreasing in j and counted 
with multiplicity), one has 

2 

(28) UJj ~ C m)k ( jfjjjq ) as j -> oo, 

where C m ,k depends on the dimensions m and k only. 

Considering the Jacobi operator 3 for a minimal submanifold K, it is easy to see from ()26f) 
that, since 3 differs from — A^ only by a bounded quantity, we have the same asymptotic formula 
for its eigenvalues ((J,i)i, and thereby 

(29) m ~ C m , fc (y^j) " as I - oo. 

In the following, we let (^i)i (resp. (<fj)j, (tpi)i) denote a base of eigenfunctions of — A% (resp. 
of — A^, 3); normalized in L 2 (K) (resp. in L 2 (K; NK)), namely the set of functions (resp. 
normal sections of K) satisfying 

-A K (f>i = pifa; —Aftipj = uijipj] Ztpi = Viipi, = 1,2, 

Finally, using the eigenvalues (pj)j and (pi)i, one can express the L 2 norms, or the Sobolev 
norms of linear combinations of the <^'s and the ipfs. In particular, if / = ^2jOij(pj, and if 
g = YliPl'&l are an L 2 function and an L 2 normal section of K, and if L\ = ]C Q c Q (y)<9^ , 
L2 = ^2 a c a {y)C^^) a are differential operators of order d with smooth coefficients acting on 
functions and normal sections respectively, then one has 



(30) || W||£ 2(K) < C Ll £(1 + ppa]; II^H^atk) < C L , £(1 + | W | d )/3; 



An estimate similar to the latter one in (|30p holds by replacing the /U/'s by the Wj's, namely if 
9' = ZjPjVj, then \\L 2 g>\\ 2 LHK;NK) < C L2 £,-(1 + kf)($) 2 - 



3 Approximate solutions to ( |P £ | ) 



In this section, given any positive integer /, we construct functions uj £ which solve (P e ) up to 



an error of order e . We will find approximate solutions of ( P s ) in the following form 



(31) Xe(\<\)( w {(' + *(ey), Cn+i) +e^i(ey, C'+$N), Cn+i) + - • •+e J u;/(ey, C'+$N), Cn 



+1/ 
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where 3>(ey) = &o(ey) + • • • + e 1 2 $/_2(ey) and where the cutoff function \e satisfies the 
properties 

( Xe(t) = l for t € [0, 

(32) J Xe(*) =0 far t € [fe-V-7], 

I |xS°(*)| <a^, i en. 

Here 3>o> • • • ! are smooth vector fields from K into NK, while wx,...,wj are suitable 
functions determined recursively by an iteration procedure. For doing this we choose a system 
of coordinates in a neighborhood of <9£\ for which the new metric coefficients can be expanded 



in powers of e, see Lemma 13.21 below. In this way we can also expand ( P e ) formally in powers 
of e and solve it term by term. The functions (lOj)i will be obtained as solutions of an equation 
arising from the linearization of ([3]) at Wq, while the normal sections will be determined 

using the invertibility of the Jacobi operator. Notice that, by the translation invariance of 
([3]), the linearized operator possesses a non-trivial kernel, which turns out to be spanned by 
{d^WQ, . . . , d^ n wo}. The role of <&o, . . . is to obtain at every step orthogonality to this 
kernel and to solve the equation using Fredholm's alternative. 

The method here is similar in spirit to the one used in [38] except for the fact that, working 
in higher dimensions and codimensions, more geometric tools are needed. Therefore, we will 
mainly focus on the new and geometric aspects of the construction, omitting some details about 
the rigorous estimates on the error terms, which can be handled as in |38j . 

3.1 Choice of coordinates near dfl e and properties of approximate solutions 

Let T : U -> dU, where U = Ui x U 2 Q R k x R n is a neighborhood of in R N ' 1 , be a 
parametrization of <90 near some point q G K through the Fermi coordinates (y, Q) described 
before. 

Let 7 G (0, 1) be a small number which, we recall, is allowed to assume smaller and smaller 
values throughout the paper. Then for e > we set 

B £n = {x G R n + +1 : |a?| < e" 7 } . 

Next we introduce a parametrization of a neighborhood (in Q e ) of | G dtt e though the map T e 
given by 

(33) T £ (y,C',Cn+i) = ^T (ey,eC) + Cn+i^ey, e C), x = (y,C'Xn+l) G ~Ul x B e,j, 

where ey = y and where v{ey,eQ') is the inner unit normal to dSl at To(ey,eC')- We have 

dT £ dT , . dv , dT £ dT , A . cV , 

^ = ^(ey,eC ) + <n^ Wa (ey,e( ); _ = — fo,,* ) + ^. +1 ^(^C )■ 

Using the equation 

(34) du x [v]=H(x)[v], 

we find 
(35) 

= [ld + eC n+1 U(ey,eO] -^(ey,e('); = [id + e( n+1 H(ey,et')] ^(ei^C) 
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Differentiating T e with respect to ( n +i we also get 



(36) 



+1 



Hence, letting g^B be the coefficients of the flat metric g = g £ (we are emphasizing the role of 
the parameter e in the entries, which is due to the dependence in e of the map T E ) of R N in the 
coordinates (y, Cn+i) ; with easy computations we deduce that 
(37) 

9af3{y,Cn+i) = g a(3 {ey)+e( n+1 (H aS g S(3 + Hpsgta) (ey)+e 2 ( n+1 2 H a5 H a pg ScT (ey), 



(v,0; 



(38) 



9aN — 0; 



9NN 



1. 



Using the parametrization in (|33p . a solution u of (P £ ) satisfies the equation 



(39) 



\/det p 



5b 5 AB \7det 5 - 9 d\ B u + « - u p = in x B ; 



e>7 



with Neumann boundary conditions on {Cn+i = 0}. Looking at the term of order e % in this 
equation, we will determine recursively the functions (lOj)j and ($j_2)j (defined in (|31|) ) for 
i = 1, . . . , I. The specific choice of the integer /, which will be determined later, will depend on 
the dimension N of Q, the dimension k of K, and the exponent p. For the moment we let it 
denote just an arbitrary integer. The main result of this section is the following one. 

Pro position 3.1 Consider the Euler functional J £ defined in and associated to problem 
(P e ) (for p < ^xjjrzf Then for any I € N there exists a function ui t£ : £l £ — > R t/ie 
following properties 



(40) 



I4K 



n/ )£ > in Q e ; 



Si 



on <9fL 



where Cj depends only on fi, i^T, p and I. Moreover in the above coordinates there holds 
' Vj m) n /i£ (y,0 < Cm^e-MPjiC), 



(41) 



Vr j V c n /i£ (y,C) < Cm^e-^PjiO, 
v£ m) V;V £ (y,0 < C m , /£ m e-lflP 7 (C), 



e 



m 



0,1,..., 



where Vy" 1 ^ (resp. V^' ) is any derivative of order m with respect to the y variables (resp. of 
order i with respect to the £ variables), where C m j is a constant depending only on Q, K, p and 
m, and where Pj(C) are suitable polynomials in £. 

In the next subsection we show how to construct the approximate solution and we give 
some general ideas for the derivation of the estimates in (I4ip . We refer to [38] for rigorous and 
detailed proofs. 

3.2 Proof of Proposition I3TT1 

This subsection is devoted to the explicit construction of uj >e . First of all we expand the Laplace- 
Beltrami operator (applied to an arbitrary function u) in Fermi coordinates, and then by means 
of this expansion we define implicity and recursively the functions (wi)i and the normal sections 



7« 
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3.2.1 Expansion of A 9e u in Fermi coordinates 

We first provide a Taylor expansion of the coefficients of the metric g = g e . From Lemma 12.11 
and formula (|37p we have immediately the following result. 

Lemma 3.2 For the (Euclidean) metric g £ in the above coordinates we have the expansions 
9ij = Sij + 2sQ n+l H %J + ^e 2 R lstj ( s (t + e 2 Cn+i 2 (ff 2 k + 0(e 3 \C\ 3 ); 



g aj = 2£( n+1 H aj + G(E 2 \t\ 2 ); 



g ab = 5 ab - 2eT b a {E i )C >i + 2eQ n+l H ab + e 2 R sabl + T c a {E s )T b c (Ei) ( s (i + e 2 ( n+1 2 (H 2 ) ab + Q(e 
g aN = 0; Snn = 1- 



Using these formulas, we are interested in expanding A 9e ii in powers of e for a function u of the 
form 



Such a function represents indeed an ansatz for each term of the sum in (|3ip . 

We recall that, when differentiating functions with respect to the variables y, £, we will mean 
that d a = d Va and di = 8q. When dealing with the scaled variables y we will write explicitly 
dy a , so that, if u is as above, we have d a u(ey, £) = edy a u(y, £). 

Lemma 3.3 Given any positive integer I and a function u : -IA\ x B £ ^ — > M of the form 
u(ey, C), we have 

A ge u = dlu + dl +l( - n+1 u + e[H«d Cn+1 u-2C n+l H lJ d 2 u] 

I 

(42) + e 2 [L 2 ,iu + L 2 , 2 u + L 2<3 u] + ^ e%u + e I+1 L I+1 u, 



where 



i=3 



L 2 ,m = dy^u - 4( n+1 H ia d CiV u; 
L 2 , 2 u = SCn+i^H^d^.u + 2( n+1 H ab r a b (E i )d i u - 2C, n+1 tr{H 2 )d Cn+1 u- 

E 2 fiU = (lliaal + -^Rihhl\ (All — -RmijlCmCld^U — -Rrniji^mdQU 

- CA^THE^u + 2QH ab T a b (E i )d (n+1 u, 

and where the Li 's are linear operators of order 1 and 2 acting on the variables y and £ whose 
coefficients are polynomials (of order at most i) in C uniformly bounded (and smooth) in y. 
The operator Lj + \ is still linear and satisfying the same properties of the Li 's, except that its 
coefficients are not polynomials in £ ; although they are bounded by polynomials in £. 
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Proof. The proof is simply based on a Taylor expansion of the metric coefficients in terms 
of the geometric properties of dO, and K, as in Lemma 13.21 Recall that the Laplace-Beltrami 
operator is given by 

A 9s = — i= d A ( ^g~ £ g £ AB d B ) , 
Vdetg £ 

where indices A and B run between 1 and N. We can write 

A, £ = 9e AB d\ B + {d A g £ AB ) d B + \ d A ( log det g £ ) g £ AB d B . 
Using the expansions of Lemma 13.31 we easily see that 

gAB q2 abU = Q2^ u + dl +lCn+ U - 2 

+£2 { 8 ha + {3(n + l 2 (H 2 )ij - I RmijlCmCl) d^u - A^H^u} + 0(e 3 \(\ 3 ). 
We can also prove 

Vdet 5e = 1 + eCn+lK + l^RmiilCmCl + \e 2 (yRmaal + T C a (E m )T a c (Ei)\ CmCl 

+ e 2 \^Cn+i 2 {Kf - ( n+ itr(H 2 ) + 2Cr+iCiH ab Ti(E i ) - C^r^r^)} 

+ 0(e 3 |C| 3 ), 

which gives 

log Vdet^ = e( n+1 HS + e 2 Ucn+xCiH^Ei) - Cn+i 2 tr(H 2 ) - CiC;^ (i^T^)} 



+k 2 R miil ( m ( l + \e 2 [R maal + T c a (E m )T a c (E l ) CmG + C(^ 3 |C| 3 )- 



Hence, we obtain 



d A (log Vd^) 9 AB d B = e 2 ^2Q n+l H ab T a b {Ei) - (^(E^iE,) + \R mhh iCi + RiaalCi\d-, 

+e^a Cn+1 n + e 2 |20^r^(ii; / )-2C n+ itr(i7 2 )|a ? , i+lU + O( e 3 |C| 



iU 
3|A|3\ 



Collecting these formulas together, we obtain the desired result. ■ 

Remarks 3.4 (a) The term of order e in the expansion of A g u in (|42p depends on the fact that 
<9f2 has an extrinsic curvature in M. . Such a term does not appear in the analogous expansion 
for the mean curvature of tubes condensing on minimal subvarieties of an abstract manifold, see 
Proposition 4-1 in \37j (where the small parameter p is the counterpart of our parameter e). 

(b) For later purposes, see for example Lemma \6.1\ it is convenient to analyze in further 
detail the operator L3 in ([42 \ , and in particular the coefficients of the second derivatives in the 
y variables. It follows from the above expansions that the coefficient of <9| ^ in L3 is given by 

2 ( Ci^a( E i) - Cn+lHab , 
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3.2.2 Construction of the approximate solution 



We show now how to construct the approximate solutions of §P £ \) via an iterative method. Given 
1 — 2 smooth vector fields <3>0; • • • > we define first the following function uj j£ on K x M n+1 , 
see (ED 



ui, £ (y, = wo(C' + $(y), Cn+i) + ewi(y, C + *(y), Cn+i) + • • • + ^wiiy, C' + $(y), Cn+i), 

where $ = $o+ e ^ > i + - • -+£"*~ 2 3>/-2- In the following, with an abuse of notation, we will consider 
uj >£ (and Wo, . . . ,wj) as functions of the variables y and £ through the change of coordinates 
V = ey- 

To define the functions (vuj)j and we expand equation (j39|) formally in powers of e 

for u = uj^ (using mostly Lemma I3.3|) and we analyze each term separately. Looking at the 
coefficient of e in the expansion we will determine w\, while looking at the coefficient of e J we 
will determine Wj and 3>j-2i for j = 2, . . . , I. In this procedure we use crucially the invertibility 
of the Jacobi operator (recall that we are assuming K to be non-degenerate) and the spectral 
properties of the linearization of ^ ai wq. 

• Step 1: Construction of w\ 

We begin by taking 1 = 1 and $ = 0. From Lemma 13.31 we get formally 

-Ag s ui !£ + ui )£ -u^ >£ = -A R n+iu> + wq - Wq + e ^-A r+ itui + w\ - pw^wi 

- e [H%d Cn+1 w - 2Q n+1 H ij d 2 ij wo] + 0{e 2 ). 

The term of order 1 (in the power expansion in e) vanishes trivially since wq solves (|3|), and in 
order to make the coefficient of e vanish, w\ must satisfy the following equation 

(43) £0^1 = H^dc_ n+1 WQ - 2( n+1 H ij df j w , 

where Co is the linearization of ([3]) at wo, namely 

J -A Wl + w x - pw p - l Wl = H%d Cn+1 w ~ 2( n+1 H tJ df jW o, in 
1 9w -t = 0, on{C„ + i=0}. 



Since Co is self-adjoint and Fredholm on H 1 (R 1 ^ +1 ), the equation is solvable if and only if the 
right-hand side is orthogonal to the kernel of Co, namely if and only if the L 2 product of the 
right-hand side with vanishes for i = l,...,n, see Proposition I4.ll below. This is clearly 

satisfied in our case since both d^ n+1 wo and dfjWo are even in while the ^ a 's are odd in 
for every i. Besides the existence of w\, from elliptic regularity estimates we can prove its 
exponential decay in Q and its smoothness in y (see for example Lemma 3.4 in [38]). Precisely, 
there exists a positive constant C\ (depending only on O, K and p) such that for any integer t 
there holds 

(44) IV^V^OI < CiC,(l + |C|) Cl e- |cl ; (y,() eKx R n+1 , 

where C\ depends only on I, p, K and f2. 
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• Step 2: Expansion at an arbitrary order 

We consider next the coefficient of e for an integer I between 2 and /, and we assume that the 
functions w\, . . . , and the vector fields <£o> ■ ■ • j have been determined by induction in 
/. The couple (vjj,<&j_ 2 ) w iU be found reasoning as for w\: in particular an equation for $/_ 2 
(solvable by the invertibility of 3) is obtained by imposing orthogonality of some expression to 
the kernel of Co, and then wj is found again with Fredholm's alternative. 

Expanding (|39|) with u = uj j£ , we easily see that (formally), in the coefficient of e , the 
function Wj appears as solution of the equation 



CqWj = Fj(yX,w ,wi, . . .,Wj_ v $ , . . . ,$/_„) ui 



n+1. 

2 J 111 ^+ ' 



( 45 ) { 9w l' n l n\ 

where is defined by 

£$u = -An + u- pw%~ ((' + $(y), ( n+ i)u, 

and where Fj is some smooth function of its arguments (which we are assuming determined by 
induction). Our next goal is to understand the role of < £/_ 2 m the orthogonality condition on 
Fj (to the kernel of £$). In order to do this, we notice that, using Lemma 13.31 for u = ui jE , the 
function <I> (precisely its derivatives in y) appears through the chain rule when we differentiate 
u with respect to the y variables. Moreover, for testing the orthogonality of the right-hand 
side in (H5j) to the kernel of £$, we have to multiply it by the functions ^ a (C' + $(y)>Cn+i)j 
i = 1, .. . ,n, so this condition will yield an equation for $ (and in particular for 3>/_ 2 ) through 
a change of variables of the form £ i— ► £ + $(y). 

Therefore, in the expansion of A g uj t£ , we focus only on the terms (of order e 1 ) containing 
either derivatives with respect to the y variables, which we collected in £2,1, or containing 
explicitly the variables £, which are listed in £2,3- In particular, none of these terms appear in 
the first line of 

Denoting the components of $ by {& 3 )j (in the basis (Ej)j of NK), there holds 

dy a (u{y, £ + $(y), Cn+i)) = dyu{y, £ + $, Cn+l) + gf-QQ^' C + $&)> Cn+l); 



Q 2 ^ (u(y, C' + *(y), Cn+l)) = 9§ a y a U(y, C + Cn+l) + 2— ajj^tifo £ + $, Cn+l 



+ ~ — — (y, C + $ (17), Cn+i) + uz- ^z- t^t- (17, C + * (y) , Cn+i ) ; 



%o ^Cj ^Cz 



<9 2 n 5^ <9 2 u 

' G7,C' + *(!/), Cn+l)) = ^ ?i%a U(y,C' + ^, Cn+l) + — U7^7-(y'C' + ^(y), Cn+l) 



<9C/%, 



<9y a <9G<9G 
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Therefore, recalling the definition of ui jE , since dy a Wo = we find that 

d 2 & dw d& d¥ d 2 w d& d 2 w 

2 Wi 



0$J 5 2 " 



• Step 3: Determining wj and ^fn for I >2 

When we look at the coefficient of e 1 in e 2 L2^ui^, the terms containing 3?/_ 2 are given by 
d 2 & dw TT d& d 2 w ( d®> d& d 2 w ., ? _ 

^Qn+lHi a —- + T^T^~ n > v 11 1 = 2 



d 2 y a Va d(j a dy a d(jd(i \ dy a dy a dQdC l 

When we project A 9e u/ )£ — ui :£ + u \ onto the kernel of namely when we multiply this 

expression by ^jp(C' + ^(l/)>Cn+l)j s = 1, ■ ■ ■ considering the terms of order e 1 involving 
2' we nave no contribution from the first line and from £2,2 in ()42|) (with u = uj, e ), as 
explained in Step 2. Also, in (|42p . the factors of for i > 3, multiplied by e I ~ 2 ^j_ 2 wm §i ve 
higher order terms. In conclusion, we only need to pay attention to £2,1 and £2,3- 

dw 
9C 

the coefficient of e 1 ^ in the following way. 

Looking for example at the first term in £ 2 L2,3 we get 



When we multiply e 2 L2^wo{C + Cn+i) by + 3>, Cn+i), s = 1, . . . , n, we can obtain 



e 2 



_2 



c 2 



/ ( -Riaal + ^i?ifc« J (AM? + $, Cn+l^oCC' + *, Cn+lR 
/ [Riaal + \Rihhl) (Cl ~ &)diWo((',( n+1 )d s Wo((',( n+1 )d( 

I [Riaal + \Rihhl) (AwoiC, (n+l)d s W {C, (n+l)d( 

JRI+ 1 V A J 

^y ^- / (flioa, + -i^ifcw J ^*00(C',Cn+l) 5 «*"0(C,Cn+l)dC- 



— £ 

i=0 

Since wq is even in it follows by symmetry that the term of order e 1 containing 3>/_ 2 in the 
last expression is given by 

(46) - C [R saal + ~R shh A 
where we have set 

(47) C = f (d lWo ) 2 . 
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From similar arguments, the third and the fourth terms in L2^wq give respectively 

(48) ^«C */_ 2J 
and 

c r*(£? a )it(^)*/_ 2 - 

The last term in L2^wq gives no contribution since the coefficient of $f_ 2 vanishes by oddness, 
so it remains to consider the second term. Integrating by parts we find 

§*mtfi*f- a ! (md <s wo% i(j md( (+dlySt 2 dl w $l 2 f al iCj w d <s w dC if 1 = 2). 

In case 1 = 2 the quantity within round brackets cancels by oddness, therefore in any case we 
only need to estimate the first one. Still by oddness in the first integral is non-zero only if, 
either i = j and m = s, or i = s and j = m, or i = m and j = s. 

In the latter case we have vanishing by the antisymmetry of the curvature tensor in the first 
two indices. Therefore the only terms left to consider are 

£?fl«iJ*f_ 2 / +1 (sd Cs w d? CiCi wod( + Y,l^sil® l i_ 2 [ Cid^wodl^wodC. 

Observe that, integrating by parts, when s/i there holds 

/ Csd^wodf f WodC = - (id Cs wod? c w dC- 

Hence, still by the antisymmetry of the curvature tensor we are left with 

-V~i? s «$~ 2 / (id Cs wodl c .w dt. 
^ 3 J ~ i 7iR™+ 1 ^ 

The last integral can be computed with a further integration by parts and is equal to — \Cq, so 
we get 

This quantity cancels exactly with the second term in (|46p and with 

dwg 
dCs 

I fhh 



When we multiply e 2 L2,iWq{( > ' + Cn+i) by ^p-iC + Cn+i), s = 1, . . . , n, the terms containing 



r$l are given by 

<9 2 <H „ r d& 



f ° ry? /_2 /" a ^/-2 2 

Jr1 +1 dy a dy a / 



which give by oddness 



8 2 <f> j - 

Cq 2 
° d 2 y a y a ' 
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Collecting the above computations, we conclude that Ff(y, £,wq,Wi, . . . ,Wj_ v $qj • ■ • } ^1-2)' 
the right-hand side of (|45p . is L 2 -orthogonal to the kernel of if and only if $/_ 2 satisfies an 
equation of the form 

Co -j^f 1 - Rsaai^j_ 2 + Y b a {E s )Y a b {E l )^ l 1 _ 2 = Gj_ 2 {yX,w ,w u . . .,w^ v <S> , *j_ 3 ), 
\ °y a Va ) 

for some expression Gf_„. This equation can indeed be solved in $f_ 2 - I n fact, observe that 
the operator acting on 3>/_ 2 in the left hand side is nothing but the Jacobi operator, which is 
invertible by the non-degeneracy condition on K. 

Having defined 3>? 2 in this way, we turn to the construction of Wj which, we recall, satisfies 
equation (|45p . Having imposed the orthogonality condition, we get again solvability and, as for 
wi, one can prove the following estimates 

(49) iV^jfoOl < C } Q(1 + |C|) C 'e-KI; (y,() G K x R n +\ 

where C\ depends only on Z, p, K and f2. 

As already mentioned, we limit ourselves to the formal construction of the functions ui e , 
omitting the details about the rigorous estimates of the error terms, which can be obtained 
reasoning as in |38j . We only mention that the number 7 has to be chosen sufficiently small to 
obtain the positivity of u/ i£ , after we multiply n/ )£ by the cutoff function Xei see (j3T|) and ([32]) . 



4 A model linear problem 

In this section we consider a model for the linearized equation at approximate solutions which, 
for p < %^ (as we are assuming until the last subsection), corresponds to J"(ui jE ). We first 
study a one-parameter family of eigenvalue problems, which include the linearization at wo of 
©. Then we turn to the model for J' e '(ui <e ), which can be studied, roughly, using separation of 
variables. 



4.1 Some spectral analysis in M™ + 

In this subsection we consider a class of eigenvalue problems, being mainly interested in the 
symmetries of the corresponding eigenfunctions. We denote points of M n+1 by (n + l)-tuples 
Cl 5 C2, • • • , Cn, Cn+i = (C'> Cn+i), and we let 

K +1 = {(Cl, C2, • • • , Cn, Cn+l) G : Cn+1 > 0} . 

For p £ fl, ^rf^ (n— I * s ^ ne cr itical exponent in M n+1 ) we consider problem ([3]) which, wc 
recall, is 

'-Au + u = u p inM™ +1 , 
< ^ = on dRl +1 , 

u>0,u€ H^R^ 1 ). 
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It is well-known, see e.g. [51] , that this problem possesses a radial solution wq(t), r 2 = Yll=i Cf> 
which satisfies the properties 



(50) 



w' (r) < 0, for every r > 



lim^oo e r r2 Wo (r) = a„ iP > 0, lim r ->. o ^) = _i ; 



where positive constant depending only on n and p. Moreover, it turns out that all the 

solutions of ([3]) coincide with wq up to a translation in the £' variables, see [23], 
Solutions of ([3]) can be found as critical points of the functional J defined by 

(51) J(u) = \ [ (|Vu| 2 + u 2 ) - -L- [ Wr 1 : UGH 

We have the following non-degeneracy result, see e.g. [50] . 

Proposition 4.1 The kernel of j" (wq) is generated by the functions . . . , ^p. More pre- 
cisely, there holds 

J (w )[w ,W } = ~{p- l)\\w \\ 2 Hl ( R n + ly 

and 

j"(w )[v,v] > C- 1 ||u||^ 1(R , l+1) , Vu E H 1 {Rl +1 ),v J_ ^o,%wo,...,5 Cn too 

for some positive constant C . In particular, we have rj < 0, a = and r > 0, where r\, o 
and t are respectively the first, the second and the third eigenvalue of J (wq). Furthermore the 
eigenvalue r/ is simple while a has multiplicity n. 

Notice that, writing the eigenvalue equation J (wq)[u] = Xu in ^T 1 (1R" +1 ), taking the scalar 
product with an arbitrary test function and integrating by parts one finds that u satisfies 



-Au + u — pwQ u = X(—Au + u) 



m 



£>Cn 



^--0 ondR" +1 . 



The goal of this subsection (the motivation will become clear in the next one) is to study a more 
general version of this eigenvalue problem, namely 



I -Au + (1 + a)u - pu% 1 u = X (-An + (1 + a)u) in 
\§£ = on<91R™ +1 , 

where a > 0. It is convenient to introduce the Hilbert space (which coincides i? 1 (R™ +1 ), but 
endowed with an equivalent norm) 

H a =\ue H\Rl +1 ) : \\u\g = [ (\Vu\ 2 + (1 + a)u 2 ) 1 , 

[ JRI+ 1 J 

with corresponding scalar product (•,•)«• We also let T a : H a — ► H a be defined by duality in 
the following way 

(53) (T a u,v)H a = / ((Vit • Vi>) + (1 + a)uv) — p I u%~ uv; u,v E H a . 

When a = 0, the operator To is nothing but j" {wq). For a > 0, the eigenfunctions of T a satisfy 
([52]) . We want to study the first three eigenvalues of T a depending on the parameter a. 



(52) 
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Proposition 4.2 Let rj a ,a a and r a denote the first three eigenvalues of T a . Then rj a ,a a and 
T a are non- decreasing in a. For every value of a, rj a is simple and there holds 

Ia > 0; lim n a = 1. 



da 



a— »+oo 



The eigenvalue a a has multiplicity n and for a small it satisfies > 0. The eigenfunction u a 
corresponding to r\ a is radial in £ and radially decreasing, while the eigenf unctions corresponding 
to a a are spanned by functions v a ^ of the form v a> i(Q = Va(\C\)w\> * = 1> • • • > n > f or some radial 
function v a (\Q\). If u a and v a are normalized so that \\u a \\ a — H^a^Ha — 1; then they depend 
smoothly on a. Moreover we have 

\V^ a (x)\ + \V®(y a ,i)(x)\ <C ie ~^, 
provided a stays in a fixed bounded set o/R. 
Before proving the proposition we state a preliminary lemma. 

Lemma 4.3 Let r denote the third eigenvalue of J (wq)- Then, for a > 0, every eigenfunction 
corresponding to an eigenvalue A < ^ of (|52p is either radial and corresponds to the least 
eigenvalue, or is a radial function times a first- order spherical harmonic (in the angular variable 



X) with zero coefficient in Q' , and correspond to the second eigenvalue. 



Proof. First of all we notice that, extending evenly across <9R™ +1 any function u G H 1 (R r ^ +1 ) 
which is a solution of (|52[) . we obtain a smooth entire solution of —Au + (1 + a)u — pw$~ 1 u = 
A (—Au + (1 + a)u). Next, we decompose u in spherical harmonics in the angular variable 9 (we 
are using only spherical harmonics which are even in Cn+i) 



2>(KI)ne(0); CeR n+1 ,0 = 7^ 

i=0 ^' 



Here Yi te is the j — th eigenfunction of — Agn (which is even in Cn+i)> namely it satisfies AgnYi e = 
\f e Yi ie , where we have denoted by \f e the i-th eigenvalue of — A^n on the space of even functions 
in Cn+i- In particular, the function Yo,e is constant on S n and correspond to Af™ = 0, while 
Af e = n has multiplicity n. The eigenfunctions corresponding to Afg are (up to a constant 
multiple) the restrictions, from R n+1 to S n , of the linear functions in 
The laplace equation in polar coordinates writes as 

A R n + lU = A r u H ^A$nU, 

where A r = + 7 Jr- Therefore, if u = Y^Zo u i(\(\) Y i,e( d ) is a solution of (J52J), then every 
radial component m satisfies the equation 



(54) 




V--V + I l + a+^-)v) -pw%~ 1 v = in 
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We also notice that, since the space of functions {v(r)Yi^ e (9)} (for a fixed i) is sent into itself by 
the Laplace operator, every Fourier component (in the angular variables) of an eigenfunction of 
([52]) is still an eigenfunction. 

We call X a ,i,j the j-th eigenvalue of ([33]) . From Proposition 14.11 it follows that Ao,i,i = 
— (p — 1) < and that Ao,ij > t for j > 2. In fact, a radial eigenfunction of J (wq) which is 
not (a multiple of) wq itself must correspond to an eigenvalue greater or equal than t, which 
is positive. On the other hand, it follows from Proposition 14.11 that Ao,2,i = 0, and also that 
Ao,2j > t > for j > 2. Finally, since Ao,i.i > r > for i > 3, we have in addition Xo,i,j > t f° r 
every i > 3 and for every j > 1. 

After these considerations, we turn to the case a > 0, for which similar arguments will apply. 
Solutions of (|54|) can be found as extrema (minima, for example) of the Rayleigh quotient 



(55) 





W? + [ 


'l + a + 




r r , p— 1 2 

-pJ R+ r n w Q v z 




{v'Y + 


'l + a + ^ 





from a standard min-max procedure. Using elementary inequalities it is easy to see that the 
above quotient is non-decreasing in a. Therefore it follows that X a ,i,j > for j > 2, that 
A a ,2,j > t > for j > 2 and that X a ,i,j — T f° r every i > 3 and for every j > 1. This concludes 
the proof. ■ 



PROOF of Proposition 14.21 The simplicity of r\ a can be proved as in [39|, Section 3, using 
spherical rearrangements and the maximum principle. The weak monotonicity in a of the 
eigenvalues can be easily shown using the Rayleigh quotient in the space H a , as for ([53]) . 

The smoothness of a i— > r] a and of a i-> m q can be deduced in the following way. Since the two 
spaces ^(M^ 1 ) and H a coincide, and since the eigenvalues of an operator do not depend on 
the choice of the (equivalent) norms, we can consider T a acting on i? 1 (M" +1 ) endowed with its 
standard norm (independent of a). Having fixed the space, we notice that the explicit expression 
of T a is given by 



(56) T a u = [-A + l] _1 (-Au + (l + a)u-pw% 



In fact, letting T a u = q G ff 1 (M™ +1 ), taking the scalar product with any v G i^ 1 (M™ +1 ) and 
using ([53]) we find 

/ [(Vq • Vv) + qv}= / [(Vn • Vv) + (1 + a)uv) - pi w^uv, 

which leads to (]56p by the arbitrarity of v. It is clear that the operator in ()56[) depends smoothly 
on a and therefore, being rj a simple, the smooth dependence on a of r/ Q and u a follows. 
We now compute the derivative of r\ a with respect to a. The function u a satisfies 



(57) 



(1 - rj a ) (-Au a + (1 + a)u a ) = pw p Q l u a in R" +1 , 
^ = ondI^ +1 . 
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Differentiating with respect to a the equation Hu^Ha = 1) we find 



(58) 



da 



l a\\a 



du a 
da 



,u c 



n+l 



On the other hand, differentiating (|57p . we obtain 
(59) 

& (" A ^ + (1 + «K) + (1 - Va) ("A (fe) + (1 + a) + u Q ) = jwT 1 * 



9_ /dUa\ _ n 



9v 



in 
on 



Multiplying (|59l) by u a , integrating by parts and using (|58|) . one gets 



(60) 



da 



(1 " Va) 



ut > 0. 



Indeed, since T a < Id Hl mn+u, every eigenvalue of T a is strictly less than 1, and in particular 
(1 — rj a ) > 0. We now consider the second eigenvalue a a . For any a > it is possible to make 
a separation of variables, finding eigenfunctions of (152 p of the form Yi jB v aj i, where Yi >e = A, 

i = 1, . . . , n, correspond to Af™. Also, from Lemma HT3l we know that for a close to (indeed, as 
long as a a < r) every eigenfunction corresponding to a a is of this form, for some i G {1, . . . , n}. 
Therefore, if we restrict ourselves to the space of functions of the form ^(ICI)ra f° r a fixed 
i E {1, . . . , n}, the first eigenvalue for (|52p becomes simple, so we can reason as before, obtaining 
smoothness in a and the strict monotonicity of a a . 

We prove next that the eigenvalue r\ a converges to 1 as a — > +oo. There holds 



inf 



f Rn +i | Vn| 2 + (1 + a)u 2 - pwQ l v? 



u€H a J R „ +1 [|Vu| 2 + (1 + a)u 2 ] 

Fixing any 5 > 0, it is sufficient to notice that 

\Vu\ 2 + ((1 + a) - pw^' 1 ) u 2 > (1 - 5) [\Vu\ 2 + (1 + a)u 2 ] 



for every u, 



provided a is sufficiently large. This concludes the proof of the claim. 

The decay on u a , v a i and their derivatives is standard and can be shown as in 
do not give details here. ■ 



so we 



Remark 4.4 Proposition ^. S\ implies in particular that there is a unique a > such that rj^ = 0. 
Moreover, we have also 

n = Coujo; Vq = Cod h w Q , 

for some positive constants Co and Cq. 



We also need to introduce a variant of the eigenvalue problem (|52p . for which we impose vanishing 
of the eigenfunctions outside a certain set. For e > and for 7 e (0, 1) we define 



(61) 



D 



e>7 



{x £ 



pn+1 



x < e 
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and let 

Hi = {u G H l (B £n ) : u{x) = for \x\ = e~ 7 } . 
We let H a:£ denote the space H £ endowed with the norm 

\Hl,e= [ [\Vu\ 2 + (l + a)u 2 ]; ueHl 

and the corresponding scalar product (•, •)<*,£• Similarly, we define T afi by 

(T U)S u, v) a>s = / (Vu • Vv) + (f + a)uv - pv%~ uv 



u,v G H a , E . 

The operator T Qj£ satisfies properties analogous to T a . We list them in the next Proposition, 
which also gives a comparison between the first eigenvalues and eigenfunctions of T a and T a ^ £ . 

Proposition 4.5 There exists £q > such that for e G (0, £o) the following properties hold 
true. Let r) aE , o~ ay£ and t Qj£ denote the first three eigenvalues ofT aE . Then r/ atE , cr ay£ and T a ^ £ 
are non- decreasing in a. For every value of a, rj ajE is simple and —^f- > 0. For a sufficiently 
small, G a)£ has multiplicity n and > 0. The eigenfunction corresponding to r/ Qj£ is 

radial in ( and radially decreasing, while the eigenfunctions corresponding to <7 aj£ are spanned by 
functions v aj£: i of the form v a:£j i(C) = v a ,e(\C\) m > * = 1> ••• > n > f or some radial function v a ,e(\C\)- 
The eigenvector u a , £ (resp. v a>E j), normalized with ||u a)£ ||ff = 1 (resp. ||va,e,i||if„ e t = lj 
corresponding to rj aE (resp. a ai£ for a small) depend smoothly on a. Moreover for some fixed 
C > there holds 

(62) |V(V, e (C)l + |V(Ve,i(OI <Cie~$, fori = 0,...,n; 



(63) \J] a — T] a ,e\ + \\u a — Ua,e\\jjlQK n + ls j 4" \o~a ~ Ca,e| "I" \\ v a,i ~ ^a,£,i||j^i(ig>™+ 1 ) ^ Cc c 5 

provided a stays in a fixed bounded set o/R. TTte functions u a)£ and v aEs i in this formula have 
been set identically outside B E ~. Furthermore, r at£ > r a > r for every value of a and e. 

The proof follows that of Proposition 2.3 in [3D], and hence we omit it here. It is still based on 
some elementary inequalities and on the Rayleigh quotient. The quantitative estimates in f|63|) 
can be deduced using cutoff functions and the Green's representation formula for the operator 
-A + (l + a) in 

As a consequence of this proposition (taking a = 0) we obtain that, if (for e small) u G H £ 
has no Fourier components (in 8) with indices less or equal to n, then (To >E u, u)o )E > ^(u,u)o >e . 
Equivalently, there holds 
(64) 

P <~\\(\y<(l-~) (-Au + u)ud( forany U = £ Uj (\(\)Y i>e (e), u G 

J B a JB E>1 i=n+l 
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4.2 A model for J'J(u Ij£ ) 

In this subsection, using the analysis of the previous one, we construct a model operator which, 
up to some extent, mimics the properties of J' e '(ui >E ), and for which we can give an explicit 
description of the spectrum. Although the related construction in [38J is a particular case of the 
one made here, the general spirit is quite different, and is more geometric in nature. 

First of all, we choose an orthonormal frame (Ei)i as before, and we define a metric g on NK 
as follows. For v £ NK, a tangent vector V £ T V NK can be identified with the velocity of a 
curve v(t) in NK which is equal to v at time t = 0. The same holds true for another tangent 
vector W S T V NK. Then the metric g on NK is defined on the couple (V, W) in the following 
way (see [21], pag. 79) 



g(V,W) = g(^V,^W) + 



D N v, D N w, 



\t=Q, T, — \t=0 



dt dt i N 



In this formula n denotes the natural projection from NK onto K, and denotes the (normal) 
covariant derivative of the vector field v(t) along the curve tt v(t). In the notation of Subsection 
Owe have that, if v(t) = v j (t)Ej(t), then 

D N v dvHt) , . , M „,/ dv(t)\ „ 

Therefore, if we choose a system of coordinates y on K and then a system of coordinates on 
NK defined by 

(y,()£R k xR n \ > ( j Ej(y), 

we get that 

9ai(y, = gM + VV (v&Ek, v%ej) n = u(y) + tVti (%) /?] (%) , 

and 

9ai(yX) = C J f3 l 3 (da); %(y,C) = %, 

where we have set = We notice also that the following co-area type formula holds, for 
any smooth compactly supported function / : NK — > M 

(65) / fdVg = [ ( [ f(£)d() dVj{y). 

J NK JK \JNyK ) 

This follows immediately from the fact that det<? = detg, which in turn can be verified by 
expressing g as a product of three matrices like 

Id \( g \ ( Id 
Id J \ Id J \ CP Id 

the first and the third having determinant equal to 1. 
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Having denned the metric g, we express the Laplacian of a function u defined on NK with 
respect to this metric. In Fermi coordinates centered at some point q £ K, using (|16p . (|17p and 
(jT2D , it turns out that (for y = 0) 

(66) Ag U = d^u + d^u. 
Next we define the set S £ as 

S £ = {(v,(n+l) G ^ x M+ : (M 2 + Cn+i 2 )^ < £~ 7 } , R+ = {Cn+1 : Cn+1 > 0} , 

where NK £ stands for the normal bundle of K £ (in f2 e ). We next endow S* e with a natural 
metric, inherited by g through a scaling. If R £ denotes the dilation x i— > ex in (extended 
naturally to its subsets), we define a metric y e on S £ by 

In particular, choosing coordinates (y,C) on iVJ^ via the scaling (y, C) = £(y,C)> one easily 
checks that the components of g £ are given by 

(9e)ab(y,v) = (g)^(sy) + s 2 vVpl (d*) (ey)/?j (%) (ey), 

(9e)ai(y,v) = (<9a) (ey); (g £ )ij(y,v) = 

and also 

Therefore, if u is a smooth function in S £ , it follows that in the above coordinates (y , (' ', Cn+i) 
(at y = 0) 

(67) A^ u = a a 2 aU + ^ + a| n+iCn+i n. 

In the following, to emphasize a slow dependence of a function -u in the variables y, we will 
often write u(y,() = u(ey,C) (where, we recall, £ = (C',Cn+i))> identifying with an abuse of 
notation the variable y parameterizing K £ with y, parameterizing K. In this case we have that 
(at the origin of the Fermi coordinates) 

(68) A §E u = e 2 d^u + dlu + d} n+lU+1 u. 

For later purposes, we evaluate Ag s on functions with a special structure. In particular, if 
we deal with a function u of the form u(y, = <fi(y)v(\C\), we have that 

(69) A- ge u = e 2 (A K 0(y))v(\(\) + <P(y)\v, 

and if instead u(y, C) = v(\C\)tp h (y)j^ for some smooth normal section tp = tp^E^, then we find 

(70) A~ ge u = e 2 (A N K ^)\y)^v{\C\) + ^(y)A c («(|C|)|[) • 
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Now we introduce the function space Hs e denned as the family of functions in H l (S £ ) which 
vanish on {\v\ 2 + Cn+i 2 = £~ 2 " 1 }, endowed with the scalar product 

(71) (u, v) Hse = [ (V §E n • V Se v + uv) dV~ 9e . 
We consider next the operator Tg e : Hs e — ► H$ e defined by duality as 

(72) (T Se u,v)h Se = J {y~g e u-V~g e v + uv-pwl~ l {\Q\)uv}dVg e , 

for arbitrary u,v G H$ E ■ Our goal is to characterize some of the eigenvalues of T$ E , with the 
corresponding eigenfunctions. 

For simplicity, if u a ^ £ , v a , £ ,i, Va,e and ct Qj£ are given by Proposition 14,51 recalling our notation 
from Subsection 12.21 we also set 



(73) u ji£ = u £ 2 pjt£ ; vi^i = v e 2 Wii£ii ; rj j>£ = rj £ 2 p . jE ; a ij£ = <r £ 2 Uu£ . 

We also assume that these functions are normalized so that 

K-lVe = Ib„ s (\^3,e? + (1 +£ 2 P>£ £ ) = 1; 

HaWl^s = Ib,, s (|V^, £> i| 2 + (1 + e^vf^ = I. 
After these preliminaries, we can state our result. 



(74) 



Proposition 4.6 Let eq,e be as in Proposition \4-5\ Let A < ? be an eigenvalue ofTs £ . Then 
either A = r/j j£ for some j, or A = a^ £ for some index I. The corresponding eigenfunctions u are 
of the form 

(75) u( y ,c)= Yl ajfij^yhjAO + Y A¥>K £ 2/He,i(C), 

where (y,() denote the above coordinates on S £ , and where (aj)j, are arbitrary constants. 
Viceversa, every function of the form (|75|) is an eigenfunction of T$ £ with eigenvalue A. In 
particular the eigenvalues of T$ £ which are smaller than j coincide with the numbers (r]j !£ )j or 
(o~i t£ )i which are smaller than j. 

Proof. The proof is based on separation of variables and the spectral analysis of Proposition 
14.51 Integrating by parts, one can check that the eigenfunction u of T$ £ satisfies the following 
equation 



(76) 



(1 - A) (-Ag e u + u) -pwl l {C)u = in S £ , 
8^=0 on { Cn+1 = 0}. 



As before, we can extend u evenly in Cn+ij to obtain a smooth solution of the differential equation 
in ([7SD in the set {(v, ( n+1 ) G NK £ x E : (|t;| 2 + Cn+i 2 )^ < £~ 7 }- Hence, fixing y G K £ , we can 
use Fourier decomposition in the angular variable of £, and we can write 



(v,C)=X>(y,|Cl)M0)> 



1=0 
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where 6 = X £ S n , and where }/ )e is the Z-th spherical harmonic function which is even in Cn+i- 
We now decompose u further in a convenient way as 

(77) u = u + u 1 + u 2 , 

where 

u = -^=u (y,\C\); ui= E «i(»,ICI)neW; «a= E «*(v,ICI)n.W- 

vP I l=l,...,n l>n+l 

Integrating by parts, the last formula, together with ([65]) . ([69]) and (fTPj) (recall that Yj e for 
Z = 1, . . . ,n are linear combinations of j& on 5 n , /i = 1, . . . , n) easily imply that (u.i,Uj)H Se = 
for 2 7^ j and that (Ts s u i ,y L j)H Se = for i ^ j, namely that T$ £ diagonalizes with respect to the 
above decomposition ([771) . 

We begin by considering the action of T$ E on u . Using a Fourier decomposition of u (y, |C|) 
through the eigenfunctions {<pj)j of the Laplace-Beltrami operator on K £ we set 

oo 

u (y,\(\)=J2<t>i(zy)uj(\(\)- 

3=0 

By (|69p we get immediately that for any j 

A §e (^N)%(|CD) = (e 2 A^ + A c )(^( £ y)^(|C|)) = (A c - e 2 p^ey)^)- 

As a consequence we find that u G JET* satisfies the following partial differential equation in 
-B £)7 , with Neumann boundary conditions on {Cn+i = 0} 

oo 

-A^Mo+Mo-^r'dciK = E^-(^) (- A c^'(ici) + (i+£ 2 Pi)%-(ici) -^r'dciK-dci) 

3=0 

From this formula it follows that if T$ e u = Xu for some A, then by the orthogonality to u l5 u 2 
we have also Ts e u = Au , and each of the components Uj (which are radial in £) satisfies the 
eigenvalue equation T £ 2 p . £ Uj = Xuj in H £ 2 p . £ with the same value of A, where we are using the 
notation of Subsection 14,11 Using the same terminology, we can further decompose uj as 

^i(KI) = a j u j,e +^',£ with a>j S K and with (uj t£ ,Ttj }£ ) £ 2 p . £ = 0. 

From the spectral analysis carried out in the previous subsection it follows that if A < ? (and e 
is sufficiently small), then Uj E = for every j, and A = i]j )£ for some set of indices j. 

We now turn to the evaluation of Tg e on u^. Similarly as before, expanding with respect to the 
eigenfunctions of the normal Laplacian we can decompose u x in the following way 

- Ci 
Mi(y,C) = EE^'CIVmN)^ 

l>0 i=l <^< 
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and from ([70]) we deduce that 



A 9e ( E^lKDHiNil 



vi=l 



ICI 



(A C - £ %)(^^(IC|)^N)|| 



vi=l 



As a consequence we find that also 



Z>0 i=l 



-Ac ( Wl)|r) + (1 + b%)^(ICI)4 -K _1 (IC|)^(IC|)^ 



Hence, by the spectral analysis of the previous subsection, reasoning as for u we deduce that 
if n x satisfies Ts e Ui = Au x with A < |, then uj(|C|)in = u i,e,i) an d hence it follows that A = <7; )E 
for some set of indices I. 

Finally, we turn to u 2 . Proceeding as for the definition of the metric g (and using the same 
notation), we can introduce a bilinear form q (semi-positive definite) on T NK defined by 



s(V,w) 



D N v , D N w 
t=o 



dt 



dt '* °/ A 



Using again a scaling in e, we can also introduce the following bilinear form on S e 

q e = \(R £ )*Q <S> d( n+ i 2 . 
The components of this form in the above coordinates (y, £) are given by 

(Qe)ab(y, v) = eWpi (da) {ey)$ (%) (ey); (g £ ) ai (y, v) = ev^) fa) (ey); 

(3e)ij(y,v) = Sij] (q £ )nN = 1; (fle)iVa = 0. 

We then define by duality the operator T £ through the formula 

z (Vg e u,V~ gs u)+u 2 -pwl-\\C\)u 2 



(1 e u,u) Hs£ ■-- 



dV Se . 



Moreover, computing the pointwise action of T e integrating by parts, reasoning as for the deriva- 
tion of ([68]) . and using (f65j) . one finds that 

(78) (% £ u,u) HSe = j j (-uA ( u + u 2 -pwodtflP- 1 ) d( dVg e (y), u G H Se , 

where we have set <? £ = ^(R E )»g and S VjE = |(u,Cn+i) G N y K E x R + : (\v\ 2 + Cn+i 2 ) 2 <£~ 7 }- 
Hence, using ([65]) (with the scaled metric g £ ), with u = u 2 an d (ESI) we find 

p [ wl~ l ^dv~ ge =p [ [ [ w^y 2 ) dV- 9e (y) < (i - J) / / (-^Acm 2 + m1) 
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Since r < 1 (being an eigenvalue of J"(wq) < Id H i^n+i^), we deduce that 

(T Se u,u) Hse = ( c Z £ u,u)h Se + I [{g £ - q £ )(V h u,V - 9e u) + u 2 ] dVg e 

^ W [fle(Vg £ n, V 5e «) + u 2 ] dV~ E + / [{g £ - Q £ )(V gE u, V h u) + u 2 ] dV~ 9c 

> 2 Nik- 

If follows that there are no eigenvectors of the form u 2 corresponding to eigenvalues smaller than 
~. This concludes the proof. ■ 

Remark 4.7 For later purposes, it is convenient to consider a splitting of the functions in Hs £ 
which is slightly different from the one in <\77\i . If Uq, U\ and u 2 are as above, with 

uo = 1 %2 ( i ) j( £ y)U'j(\(\y, Mi =5^5^«!(ICI)pm( e iO|7|> 

j>0 l>0 i=l 

for some real sequences (ttj)j, ((3i)i, we can write 

%'(ICI) = Oj%>(ICI) + %e(ICI), with (u j:£ ,u ji£ ) £ 2 p . )S = 0; 



^(ICI)|7| = PlVl,e,i(C) +^,e(ICI)|7| := A«Z,e,i(C) + Vj,e,i(0> w ^ Ke,i, U^,*)^^ = 0. 



||r = PlVl,e,i{() +Vj,e(|CI)j|j 

./Vow we se£ u = Uq + Ui + u 2 , where 



u = ^2aju j)£ (\c\)(f>j(£y)) ui = ^A^, e ,i(C)yK £ y); 



J'=0 



TTien &j/ ([71]) one can c/iecfc f/iai (ui,Uj)# ge = /or i 7^ j, and i/iai 

(79) \\u\\ 2 Hse = \\uo\\ 2 HSe + \\ Ul \\ 2 Hse + \\u 2 \\ 2 H3e = - + - ^/f + ||u 2 

j=o 1=0 
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(80) (T Se u,u) Hse = Y,Vj,ea 2 + Y,°i,eP? + (Ts £ u 2 ,u 2 )H Ss ; (T Se u 2 ,u 2 ) HsB > CKH^, 

j=0 1=0 

for some fixed positive constant C . 

From the last proposition we deduce the following corollary, regarding the Morse index of Tg e 
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Corollary 4.8 Let 7 E (0, 1), and let Ts £ : Hs e — ► -ffs e &e defined as before. Then, as e tends 
to zero, the Morse index ofTg e satisfies the estimate 



where a is the unique real number for which rja = (see Remark \4-4\ )- 

Proof. From Proposition 14.61 we have that the Morse index of T$ e is equal to the number 
of negative rjj^s. By the estimate in (I63p . this number is asymptotic to the number of j's for 
which Tj e 2 p . is negative. Therefore it is sufficient to count the number of eigenvalues pj for which 

k 

e 2 pj is less than a. By the Weyl's asymptotic formula, see [32], we have that pj ~ C/. ( V J( K ^ ^ 2 
so the conclusion follows immediately. ■ 



k 

a \ 2 



5 Accurate analysis of the linearized operator 

In this section we first compare J £ (uj t£ ) to the model operator introduced in the previous one. 
A naive direct comparison will give errors of order e, see Lemma 15.11 and Corollary 15.31 but 
sometimes we will need estimates of order e 2 . Therefore we will expand at a higher order the 
eigenvalues (of the linearized operator at ui t£ ) close to zero with the corresponding eigenfunc- 
tions, to get sufficient control on the errors. Finally, using these expansions, we will define 
a suitable decomposition of the functional space for which the linearized operator is almost 
diagonal. 

5.1 Comparison of J'J(ui j£ ) and T$ E 

We define first a bijection T £ from S £ into a neighborhood of K £ in Vt £ in the following way. 
Given the section = $0 + £< ^i + ■ ■ ■ + £ I_2 ^i-2 in NK constructed in Section [3l for any 
{y, Cn+i) 6 S £ , v E N y K £ , Cn+i e R+, we set 

f E (v,(n+i) = exp^(v + ^(sy)) + C n+1 u(exp^(v + ^sy))) . 

Then we define the set S £ C fi e to be 

S £ = T e (S e ), 

endowed with the standard Euclidean metric induced from WL N . For u E Hg s , we define the 
function u : T, £ — > R by 

u(z) = u (f ~ l {z)j , z E S e , 

and letting A e to be the map it 1— > u, we define 

H Ss =A e (H Ss ). 

H% £ has a natural structure of Hilbert (Sobolev) space inherited by i!f 1 (J7 e ), and we denote by 
('i')h-e j II ' 1 1 He the corresponding scalar product and norm. More precisely, we can identify 
the space H-£ e with the family of functions in which vanish identically in ft £ \ E e . 
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We introduce next the operator T^ £ : — > Hy, s denned as the restriction to Hj] E of J £ (ui >£ ) 
which, using the duality in Hy, s , has the following expression 

(81) (Tx s u,v)hv 5 = (Vu-Vv + uv)-p u p I ~ 1 uv = (u,v) Hj , s -p u^uv. 

Fixing these notations and definitions, following the arguments at the beginning of Section 
4 in [38 j one can easily prove the following result. 

Lemma 5.1 Identifying the functions in Hs^ with the corresponding ones in H% e via the map 
A £ , for e sufficiently small one has 

(u,v) Hse = (u,v) Hse +0(£ 1 ~ 7 )IM|// s JM|// Se ; 

(Tx e u,v) Hse = (Ts E u,v) Hse +0(E 1 " y )\\u\\ Hs J\v\\ Hse . 
with error 0(e 1-7 ) independent ofu,v G i?s e - 

We introduced the operator T^, £ because it represents an accurate model for J £ (ui i£ ). In 
fact, since most of the functions we consider have an exponential decay away from K £ , it is 
reasonable to expect that the spectrum of J £ (ui )£ ) will be affected only by negligible quantities 
if we work in Hs e instead of JET (fl e ). More precisely, one has the following result (we recall the 
definition of r from the previous section). 

Lemma 5.2 There exists a fixed constant C , depending on £1, K andp such that the eigenvalues 
of J'J(ui )E ) and Ts e satisfy 

IXjiJ'Jim^)) - A^TsJl < Ce'c^y, provided X^j'^u^)) < T -. 

Here we are indexing the eigenvalues in non- decreasing order, counted with multiplicity. 

We omit the proof of this result because it is very similar in spirit to that of Lemma 5.5 in [39] . 
This is based on the fact that the number of the eigenvalues of T$ e which are less or equal than 
fr is bounded by e~ D for some D > (see Proposition 14.61 and the Weyl's asymptotic formulas 
in Subsection I2.2p . together with the exponential decay of the eigenfunctions of J"(ui j£ ), which 
can be shown as in [39], Lemma 5.1. 

As a consequence of Lemmas 15.11 and 15.21 we obtain the following result. 

Corollary 5.3 In the above notation, for e small one has that 

(82) IW'e'M) - X J( T S,)\ < Ce l ~\ provided A^J^u/,,)) < T -. 

Using Proposition 14 . 61 and Corollary 15.31 it is possible to obtain some qualitative information 
about the spectrum of the linearized operator J £ {uj £ ). However, this kind of estimate is not 
sufficiently precise by the following considerations. First of all, since the eigenvalues of T$ s 
can approach zero at a rate min{e 2 ,e fc }, the estimate (|82j) need to be improved if we want to 
guarantee the invertibility of j"(uj E ). Furthermore, it would be natural to expect that the 
Jacobi operator (and its invertibility) plays some role in the expansion of the eigenvalues, and 
this is not apparent here. 

On the other hand, Lemma 15.21 gives an accurate estimate on the eigenvalues of J £ (ui i£ ) in 
terms of those of Ts e , so it will be convenient to analyze Ts e directly. 
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5.2 Approximate eigenfunctions of T^ e 

In this subsection we construct approximate eigenfunctions to the linearized operator at the 
approximate solutions u/ >£ . By the reasons explained at the end of the previous subsection, 
we need a refined expansion of the small eigenvalues of Ts e , and in particular here we want to 
understand how the <7; )£ 's change when we pass from Ts s to Ts £ . 

It is sufficient here to take 1 = 2, because the terms of order higher than e 2 do not affect 
the expansions below. As for the construction of the approximate solutions uj ie , we proceed by 
expanding the eigenvalue equation formally in powers of e. By the construction of it2, £ , formally 
the following equation holds 

-Ag E U 2 , £ + U2,e ~ U P 2 e = 0{e 3 ). 

Using Fermi coordinates as in Section [3] and differentiating with respect to C,h-> we get 
(83) - d h (A 9i: U2,e) + d h u 2 , E - pvP^d^e = 0(e 3 ). 



From the general expression of the Laplace-Beltrami operator, see formula (|14p . we can easily 
see that 

d h (A gE u) = A ge (8 h u) + d h gf B d AB u + d h {d A gf B )d B u 
(84) + ^ B d 2 hA (log(detg £ ))d BU + ^d A (log(detg e ))(d h g^ B )d B u. 

Let us now consider the second term on the right-hand side of (|84|) : dividing the indices this is 
equivalent to 

d h g^d 2 ljU + 2d h gfdlu + d h gfd ab u + 2d h g* N d A d Cn+1 u. 

From Lemma 13.21 and using the fact that we get an e factor each time we differentiate u with 
respect to y a , yj, . . . , we find that 
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Similarly we get 



d h g* B d\ B u = --e 2 R lhtj Ctdf jU + 0(e 3 ) 



d h d A g^ B d B u = -s 2 R hii jdjU + 0(e 3 ); 



\g* B d 2 hA (log(detfc)) d B u = e 2 (±R illh + R iaah - r b a (Ei)T$(E h )j d iU 

+ 2H ab T b a (E h )d (n+1 u + 0(e 3 ), 

and 

\d A (log(det &)) (d h g? B )d B u = 0(e 3 ). 
Putting together all these terms we deduce that 

(85) d h (A ge u) = A gE (d h u)-^e 2 R ihtj (td ijU +e 2 {^R inh + R iaah - r b a (Ei)T a b (E h )^ d lU +0(e 3 
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To construct the approximate eigenfunctions d £ and the approximate eigenvalues n, we make an 
ansatz of the type 

t> e = (^ h (y)d h u 2 , £ (y, C + Cn+i) + e 2 z 2 (y, q) + 0(e 3 ); » = e 2 JI + 0(e 3 ), 

where the normal section ip = (ip h )h, the function z 2 and the real number JI have to be deter- 
mined. 

We notice that the eigenvalue equation J £ '(u 2i£ )v = Aw m H 1 ^^, with an integration by parts 
becomes 

- A ge v + v-p (n 2i£ ) p_1 v = A (-A ge v + v) , 

see also the derivation of ([57]) . 

For v = t> £ and A = fi, we have the following expansion 



- A 9e (^ h (y)d h u 2 , £ + e 2 z 2 (y,() 
+ ^ h {y)d h u 2 , e + e 2 z 2 (y, () - p (u^f 1 (^ h (y)d h u 2>£ + e 2 z 2 (y, () 



-A ge (i> h (y)d h u 2:£ + e 2 z 2 (y,0) + (i> h (y)d h u 2 , £ + e 2 z 2 (y,()) 



e 2 [i 



il, h {y) {-A ge d h w Q + d h w )\ + 0(e 3 ) 

= e^p^mw^dhWo + Oie 3 ). 
From (|85p we can expand the Laplacian in the last formula as 

-A ge U h (y)d h u 2>E ) = -e 2 d 2 i, h d h w Q -2e 2 d a i ) h d 2 h w dy a ^i-^ h A g£ {d h u 2 , 



+ 4e 2 Cn+iH aj dyJ h d 2 h w + 0(e 3 ) 

= -e 2 dl v J h d h w - 2e 2 d a ip h d 2 h w dy a & - i> h d h {A gs u 2;E ) 

2 

+ Ae 2 C, n+ iH aj dy a ^ h d 2 h WQ + -e 2 %l) h R iht jC,tdijWo 

- e 2 ^ h (^Riuh + R iaah - T b a {Ei)T a b {E h )^ d iWo + 0(e 3 ). 

Using (|83p jointly with the last equality, and recalling our previous notation (from Section [3]) 

Cq,u = -Au + u- pwl' 1 ^' + <J>(ey),Cn+i), 
we obtain the following condition on z 2 

ha2 



(86) 



£$>z 2 = d§ a¥a ip n d h w + 2dy a ijj h d 2 h w dy a % - -^ h RihtjCt%wo 

+ ^ (^Rillh + Riaah ~ T b a {Ei)Y a h {E h )\ diW + pji^w^dhwo 

- 2H ab F b a (E h )d Cn+1 w - 4( n+1 H aj dyJ m d 2 m w + 0(e). 
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In order to get solvability of this equation (in z 2 ), we need to impose that the right-hand side 
is orthogonal to the kernel of namely that, multiplying it by 8 s wq and integrating in £, 
s = 1, . . . ,n, we must get zero. If we do this, reasoning as at the end of Subsection 13.2.1} we 
obtain the following condition on ip 



C 3^ = Ci(iif), 



where 



Ci=p 



and where Co is given in ([I7]) . With the choices 

_ Co 



1p = 1pl 



where m is an eigenvalues of 3 with eigenfunction tpi, the right-hand side of (|86p is perpendicular 
to the kernel of and we get solvability in z^- Using the eigenvalue equation for ipi, ()86p can 
be simplified as 

+ ^RijjhdiWo - RihtjCtdfjWo - 3H ab T b a (E h )d in+1 w 



Next, we set 



-, i Co p-x 



^(y,C) = 2C- 1 \(dy a & - 2Q n+l H a3 ) d] h w 



9%(yX) = ^ [[Riuhdiwo - RihtjCtdfjWo - 3H ab T b a (E h )d Cn+1 wo)\ + d h w 2 (y,(' + Hv),Cn+i] 



and 



g$(y,0 = d hWl (yX' + $(y),Cn+i). 



We notice that, by the definitions of Co, C±, the computations in Subsection l3.2.2l and by oddness, 
the arguments of in the definitions of g\ , g\ and g\ are all perpendicular to the kernel of 
£$, and therefore go^gi and g 2 are well defined. 

Finally, with this notation, we define the approximate eigenfunction ^fi as £ times a suitable 

cut-off function of namely 

(87) 

HV,0 = Xe(\C\)\^i(V) \d h wo+eg h ^y,Q)+e 2 g^(y,Q)]+e 2 ^4{y)g h ^y 



where Xe is as in (|32|) . and, as usual, y = ey. 

A more accurate analysis, which we omit, shows that the above error terms not only are of 
order e 3 , but they decay exponentially to zero as |£| tends to infinity. Moreover, as we already 
remarked, in the above estimates one can replace with itj >e . Precisely, one can prove the 
following result. 

Lemma 5.4 If^i is given in (|87|) . then there exist a polynomial P(C) and a sequence of positive 
constants (C); ; depending on £1, K, p and I such that 



Co 



< Qe 3 P(C)e 



-ICI 
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5.3 A splitting of the functional space 



In the previous subsection we expanded in e some of the eigenvalues of Ts £ , precisely those 
which are the counterparts of the o"/ j£ 's for T$ s . Actually, Tg e possesses another type of resonant 
eigenvalues, namely the rjj^s for suitable values of j, which in principle could approach zero even 
faster. One of the differences between these two families of eigenvalues is that the eigenfunctions 
corresponding to the resonant ctz j£ 's oscillate slowly along dft £ , and this allowed us to perform 
the above expansion. On the contrary, the eigenfunctions related to the rjj^s possess only high 
Fourier modes, and therefore such an expansion is not possible anymore. Nevertheless, we can 
deal with the counterparts of these eigenvalues applying Kato's theorem, which on the other 
hand requires to characterize the corresponding eigenfunctions up to some extent. 

The purpose of the present subsection is to identify appropriate subspaces of Hs £ with respect 
to which Ts e is approximately in block form. Recalling the definitions in Proposition 14.51 i n 
(f73l) and in (f87l) (and also our convention about the range of an integer index), for 5 G (0, A;), 
C G (0, 1), we define the following subspaces 

(88) Hi = span {4>i(ey)u ij£ ((), i = 0, . . . , oo} ; 

(89) 

H 2 = span 1 = 0,.. . ,e~ 5 } ; H 2 = span (ey)v fe, J = £ ~ S + 1 ,--- ^s~ k \ ! 



(90) H 2 = H 2 ®H 2 ; H 3 = (H 1 eH 2 ) ± , 

where X 1 - denotes the orthogonal complement to the subspace X with respect to the scalar 
product in Hj} e . We have the following result, which is the counterpart of Proposition 4.2 
in [38j . The proof follows the same arguments, but for the reader's convenience we prefer to 
give details since the notation and the estimates are affected by the different dimensions and 
codimensions we are dealing with. 



Proposition 5.5 There exists a small value of the constant C > in (I89p . depending on Q, K 
and p, such that the following property holds. For e sufficiently small and choosing 5 G 
in (f89l) . every function u £ ffs e decomposes uniquely as 



u = u\ + u 2 + us, with U\ G Hi,u 2 G H 2 ,u% G H%. 
Moreover there exists a positive constant C , also depending on f2, K and p such that 

(Tx E u 3 ,u 3 ) > —^t\\u 3 \\ 2 h 
CC k 

The proof requires some preliminary Lemmas. Before stating them, we recall our convention 
about the symbol for two positive real values c and d. 

Lemma 5.6 Let u 2 = E?=^+i Pj1>? fo/)«i, e (ICI)% G H 2- Then 

Ce~ k 

Of) IHIk £ =(i + o(^))^ £ 
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Proof. By Lemma \b.l\ it is sufficient to estimate H^Hf^ • We notice that by (|26|) there holds 

-A^ipj = Z1>j + (05 - 9t)^ = W Vj + ((33 - 9\)iP)j. 
Integrating by parts, using (|7U[) and the last formula one finds that ||tt2||ff s becomes 

Ce~ k 



„ Cs~ k / n j. \ /• n 

( 92 ) + / E (E W £ y)Mici)te)-(E^ 



'MMICD 



ICI 



^1+^2, 



where 



E 



Ce" fe 



-A C + (1 + e 2 ^)) ( E Wfo)MKI) |f 



m=l 



EA^(^)% j£ (ICI)^); 



^ = £ 2 / E E ft ((» - *)^) m ( £ y)^(Ki) fer • E Atfwoj, 



(ICI) 



^2 

icir 



' j,J=e-*+l 

Looking at Al, the integral over any fiber N y K £ is non zero if and only if m = h (and by 
symmetry, when computing the integral we can assume both the indices to be 1). Then, from 
(|65p and from the orthogonality among different tp^s (which now are scaled in e), recalling that 

%,e(ICl)w = v i,e,m, At becomes 



Cc 



h E " 

j=e~ s +l 



■j ll w i,e,l|l e 2^-, e 



Ce~ r 

1 E ^ ^ l (iv^, l i 2 + (i+^)^ t i) 



Recalling the normalization ([74]) and the fact that 77^ = u>j + 0(1) (independently of j), see 
Subsection 12.21 we obtain that 



Ce~ k 



(93) 



^i = i E (l + 0(e 2 ))/3l- 



We turn now to the estimate of A2. By the orthogonality of the tpi's, using again (|65p and (|74p 
one finds __ _ 

/ U 2 2 dV~ ge = ^ Yl $H,M&(R«+i) < E ft'' 

Jbe j=e- s +l j=e~ s +l 

Working in a local system of coordinates (y, z) as in Subsection [421 it is also convenient to write 

U2 as 



Or 



«2(y,C) = E fm(y,\(\)(m, where / m (y, |C|) = E ft^fN) 



m=l 



j=e-' 5 +l 



ICI ' 
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If U is a neighborhood of some point q in K, where the coordinates y are defined, letting U £ = ^U, 
one has 



/ u 2 2 dVg e = f^ [ ( f fl(y,\C\)C!dC 



so it follows that 

Ce~ k 



(94) 



E/ / + f 2 m (yAC\)(ldc) dV- 9E {y) < [ ulav~ ge <± jh 0j. 



Now, we can write 



A 2 = e 2 [ u 2 u 2 dVg e , where u 2 = V &■((»- <K)V>,) m Mfy, e (KI)fe- 

As for u 2 , we can write u 2 = Em=i fm(z/, KDCm, where f m = Z)^e-«+l(® ~~ tymjfjiv, |C|)> and 
compute 



In conclusion, from the Holder inequality, from (|94p . covering i<C e with finitely-many W e 's we 
derive 

i i Ce~ k 

(95) |^1 2 | < e 2 ^ u|d%) 2 (jT 2 < C £ 2 ^||<B - <K|| L ~ E # 

Then the conclusion follows from f)93[) and (|95p . ■ 

In order to estimate the norm ||tt2||jf E j it is convenient to introduce an abstract result. 

Lemma 5.7 For j G {0, . . . ,£~ S }, and for a sequence {(3j)j, let us consider a function u : S e — > 
R o/ i/ie /orm 

<%>0 = E E ^(^f)(y)5m(C), 

j=0 m=l 

where y = ey, where L^y is a linear differential operator of order d with smooth coefficients in 
y, and where the functions g m (C) are also smooth and have an exponential decay at infinity. 
Then there exists a positive constant C, independent of e, 5 and (/3j)j such that 



NW £ )<^E(i + ^l>l. 



3=0 

Proof. The proof is similar in spirit to that of Lemma 15.61 but here we take advantage of 
the fact that the profile g m (C) is independent of the index j (this lemma applies in particular to 
each of the summands in the definition of *$>i, see (151 
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Using local coordinates, (|6"5"j) and the exponential decay of the g m 's, after integration in £ 
we find 



j,l=0 m,h=l 



i^- 



foi some bounded coefficients (c m h)- As for ([95]) then we find (S e ) — CllV'llifd 



i (K s ,NK s ) 
C v^e ' 



and 



the last quantity, with a change of variables and by (j30j) . can be estimated with ^j =0 (l + 
e M |/Uj| d )/3|. This concludes the proof. ■ 



Lemma 



5.8 Let U2 = u 2 + u 2 = £ PjVj(ey,Q + £ M^CICDff G #2- T/ien ; 



i=o 



choosing S G in (fSUj) , one /ias 

1 



(96) 



(l + ( £ l-7 + £ 2- f)) 



'l^0|| H l (R n+l) + E ^ 
J=0 j=e-*+l 



Proof. We first claim that the following formula holds 



(97) 

Proof of (I9TD. We write 



^ol 



j=0 



_H" 1 (IR"l +1 )' 



u 2 = n 2 ,i + n 2 , 2 := ^^^f(ey)9 m nJo(C)x £ (lCl) + E&%N>0- 

3=0 j=0 

where *$>j is the term of order e (and higher) in fyj. Reasoning as in the proof of Lemma 15.61 we 
get 

1 £ ~ S 

\\u2,i\\h Se = — ^l 32 j( l + e2 ^3+°( e2 ))\\ d niW0Xe\\ 2 HH - B n+i ) 
3=0 

= ? E^( 1 +°( £2 ~ f ))ii^oii 2 ¥1(R?+1) , 



(98) 



3=0 



where the last equality follows from the Weyl's asymptotic formula (|29p . 

On the other hand, using Lemma 15.71 the Weyl's formula and some computations, one also 
finds 

^11^2,211^ < C e 2 ^/?|(l + e 2 |^l)+C^ 4 E^M( 1 + e2 |^l) 



j=0 



j=0 



+ Ce 4 ^/? 2 (1 + \^\ + £ 2 | Mi | 3 ) < C (ei + e^+e^) 

3=0 j=0 
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By our choice of 5, the last formula reads 

(99) wil*<^E/!. 

3=0 



Finally, from (|98j) and (|99j) we also obtain 

(u2,1,U 2 ,2)h Se < — ^2 $ \ £ + ( £2 

£ i=o 

which concludes the proof of ([97]) . 

Proof of (|96fl . We write again u 2 = U2,i + ^2,2- Then, by the orthogonality relations among 
the ifij's, reasoning as in the proof of Lemma [5"Uj we get that (u 2 ,u 2 ,i)h Se becomes 

Ce~ k e~ s ,. / n j. \ , n 

e2 E E/ (E^((^-^) m ( £ ^(ici)fer ■ x £ (ici)E^( £ y)^o 

As above, with some computations we find 

Ce" fc 



(%, «2,l)fl-S e = °( e2 )ll^2||H s J|u2,l||// S£ =0(£ 2 )^ ]T ^. 



A 32 

i=0 

From Lemma 15.61 and (|99p we also find 

(u 2 ,u2,i) Hse <ci (x; (1+0(^/5?)] £ 2 -f (f: /#. 

\ 5=0 / j=0 

The result follows from the last two formulas. ■ 



Remark 5.9 From the proof of (J96J) it also follows that every function u 2 S H 2 can be written 
uniquely as u 2 = u 2 + u 2 , with u 2 £ H 2 and u 2 £ H 2 . 



Proof of Proposition 15.51 In order to prove the uniqueness of the decomposition it is 
sufficient to show that, for e small 

(100) (ui,u 2 ) Hs£ = o £ (1)\\ui\\h Ss \\u 2 \\h S£ , u\ G Hi,u 2 G H 2 , 

where o e (l) — > as e —* 0. Indeed, by Lemma 15.11 we have 

{ui,u 2 ) Hs£ = (ui,u 2 ) HSe + 0(e 1 ~^)\\u 1 \\ Hl2e \\u 2 \\ H ^, 
and since the functions dhWo, gft, g% and vi )£t i are odd in (' (and so also u 2 and U2,i), we get 

(ui,u 2 )h Se = {ui,u 2 , 2 ) Hs 
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where we have used the notation in the proof of Lemma [5 .81 Hence from the last three formulas, 
(1991) and form (1961) we deduce 



(101) («i,«2)h Bs <C(8 1 -y + e 2 - 2 k)\\u 1 \\ Hs J U2 \\ H ^, 
which implies (|100p . since 5 G (|, A;). 

To prove the second statement, it is sufficient to show that 

(102) (u3,v)h Se <^\\u3\\h s J\v\\h s /, ase^O, 

for all W3 G H3 and for all the functions v of the form 

res-* _ 

f= X] Avr( £ y)^,e,m(C)- 
(=0 

In fact, if we write u 3 = 113^ + 1*3,1 + 1*3,2 as in Remark 14.71 (with an obvious change of notation), 

00 00 
u 3 ,o = ^2,a j u jt£ {\C,\)<t>j{ey)\ u 3 ,i = ^A%,i(C)rfM> 

j=o 1=0 

from (I79p we find 

1 CO 

(103) lk3|li = ^E(«'+A 2 )+K 2 |lL £ - 

From (|79p . from Lemma |5. II and from the fact that 113 is perpendicular in Hy, s to 113^ G we 
deduce 

1 x 

— = ("3,0,U3,o)H Se = (u3,o,^3)h Se = 0(e 1-7 ) ||u 3 ||h Se [|u 3) o \\h Se < Ce 1 ' 1 ||u 3 . 

£ 1=0 

-Ce~ k 

Moreover from (|102p . choosing v = Y2i=o A ( / 9 [ n ( e 2/) u (,e,m(C)) an d using (|103p we get 



I ^ 



The last two formulas and f)103j) then imply 

/ 



(104) Ml <c- 



l>\Ce~ k 



for some fixed constant C. 

On the other hand, by (I80p we also have 
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Using the fact that crj j£ ~ (J e 2 w i ,e ~ £2 ^ k by Proposition 14.51 from (|1U4|) and the last formula it 
follows that 

(T Se u 3 ,u 3 )s £ > ^fc^r E $ + c\\ u w\\ 2 H Se ^ — =alKH?r Se - 

CC k iylCe-k + l CC k 

This yields our conclusion, hence we are reduced to prove (|102p . 
Proof of (|102p . By the form of v and by ([79]) . we have 

\Ce~ k 

(105) \\vf Hse 4E A 2 - 

1=0 

Using the I? basis of eigenfunctions of 3, we define the function ip and the coefficients 
{A}i=l,...,oo as 

17=v-fc 

2 ot oo oo 

«=0 Z=0 «=0 

so we have 

2 oe OO 

(106) iMii 2(OT) = E a 2 = Ea 2 - 

Z=0 Z=0 

Using these new coefficients /3j, we set (see ([73]) ) 

i=o i=e-*+i 
where Co is given in Remark 14.41 Hence we can write 

v -v = A 1 + A 2 + A 3 + A 4 + A 5 , 

with 

I7v- fc 

2 ot OO 

Ai= E ^r(^)ke,m(C)-%, m (()]; ^2 = E A^(ey)«o, e ,h(C); 

'=0 i = Ce- fc +l 
e~ 4 Ce- fe 

^ = -U ^ft-f i (£j/ ) C); A 4 = Y Mi( £ y)( v o,£,h-v he , h ); 

j=0 l=E~ S + l 

e~ s _ 
M = Y,Ml ( V ^,h ~ CoXe(\(\)d h WQ) , 
1=0 

and where ^Sfj is defined in the proof of Lemma 15.81 Since u 3 is orthogonal to H2, we get 
{u2,v)h^ e = 0, and so 

(107) (« 3 , v) Hj , g = (u 3 ,Ai) Hsi: + (tl 3 , A 2 )tf Ee + («3, ^3)ff Ee + (-"3, M)H-E e + («3, A 5 ) H ^ . 
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We prove now that ||^4j||jf Se is small for every i = 1, . . . , 5. From (|65p . the proof of Proposition 
6J Proposition 14.51 and (| 105j) there holds 



1*11** = ^ E ^Ik^-^illt < cc 2 (i + c 2 )ibii^ e < 1||,||^ £ , 



provided C is sufficiently small. 

To estimate A2 we can use Lemma 15.71 and some computations to find 



(108) 



_. 00 



Z=Ce- fc + l 



We now set </3 = X^=tc£- fc +i AVv- Since 3 = + ^(1), f° r an y integer m one finds 

< ((-A$)"V, ¥)l 2 (K;NK) + Cm [((~A^) m ~V, P)l 2 (K;NK) + (<P, L 2 (K ;N K)] ■ 



\Ce 



Since </? = X)z=o A^ii from (|106|) we deduce that 

2m 

~\ ^ I C \ k -2mi|,„||2 



+ o(,- 2 (™- 1 ))||^||i 2(OT) 



(109) 



< 



2m 

|\ fc e -2m + 0(£ -2(m-l) : 



E a 

«=0 



On the other hand, since in the basis (ipi)i, the function <p has non zero components only when 
I > Ce~ k , by the Weyl's asymptotic formula we have also that 



(110) 



¥)l 2 {K;NK) > < 



V°°- n m R 2 - 
2^l=Ce- k +lf J 'l Pi 



Using (|109p and the first inequality in (jllOp with m = 1 we get 



e 2 E WA 2 <(CC^+ 0£ (1)) E A 2 



2 l 



«=0 



Moreover, using (|109[) and the second inequality in (jllOp with m arbitrary one also finds 



l=Ce- k +l 



2m 



E ^ E A 2 - 



1=0 



Using (I105P , (|108l) and the last two inequalities (for the second one we take m large enough) , 
for sufficiently small C we find ||-A2||.ff s < ig|M|.ff ge - 
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Now we estimate ||-A3||i? s • Reasoning as for from (|1U5|) and (|1U6|) we get 



- 4 " 4 f ll7>H 2 

h \\ v \\Hs- 





Next, similarly to the estimate of A\, for small C we find 

\\M 2 HSe <^ c E A 2 ll«o >£ ,i-^,i||f j£ <cc 2 (i + c 2 )||,;||l, SE <^N|^. 

Finally, from Proposition 14.51 and reasoning as for A2, we obtain also 

\\A 5 \\ 2 HSe < ^Ce-e-^Y, fiil + ^OCe-H-f < Ce^e^^M^. 

1=0 

Taking (|107p into account, this concludes the proof of (|102p . provided we choose C and e 
sufficiently small. ■ 

6 Diagonalization of T^ e and applications 

In this section we study how the operator Ts e behaves with respect to the above splitting of 
H-£ £ in the three subspaces Hi^H^ and H3. We prove that its form is almost diagonal and we 
apply this analysis to study its invertibility for suitable values of e. 

6.1 Diagonalization 

Integrating by parts, we can evaluate the operator Ts e multiplying a test function by the fol- 
lowing quantity 

(111) & £ (u) = ^det g (-A g u + u- pu p f~ l 



u 

and integrating in the variables y and £ (using (|65p ). In Lemma 15.41 we studied & £ acting on 
the functions for any I fixed. In that lemma, our estimates depend on the value of the index 
I, and in general one can expect that they become worse and worse as I increases. The goal 
of this subsection is to derive estimates in terms of both e and I and, evaluating & e (u) on the 
functions 112 £ H2, we will keep track also of the terms of order e 3 and higher. 

In the following, we will sometimes omit the factor \e appearing in (|87p since this will only 
produce error terms exponentially small in e, which are negligible for our purposes. 

Lemma 6.1 There exist linear differential operators L\,L2,L^ (acting on the variables y) of 
order 1, 2 and 3 respectively, whose coefficients (independent of I) are smooth and satisfy the 
bounds 



(112) c a (L l )<C(l + \C\ c )e 
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and such that in local coordinates we have the following expression for & £ (^>i) 



6 e (*,) = e^fi^dhwo^ 

- 2e 3 (tiTliEi) - ( n+1 H ab + Cn+iH^Sab) (d§ aVb ^)d h w - e 3 ^^)^! 

(113) + eHn+lH^H^ ~ P<~ 1 ) ~ MP ~ 1K~V/^M " z 4 M( d ly a ^i)9o 

+ e 3 Li^i + e 4 L 3 ipi + e 4 mL^i + e 5 fj,iL 2 tpi, 

where Co, C\ are as in Subsection \5.2l 

Proof. As for the construction of the approximate solutions we can expand formally 
© e ( x I / ;) in powers of e and check carefully all the error terms, paying particular attention to 
the ones involving derivatives in the variables y a , which produce larger and larger terms (as I 
increases) in the Fourier modes. When we differentiate with respect to the variables the quan- 
tities appearing will be considered as coefficients (depending smoothly on with exponential 
decay) of the functions tpi or their derivatives in y. 

We recall that the functions iuo and (gi)i in (|87p are shifted in the £' variable by the (smooth) 
normal section 3>(y). Hence, when differentiating with respect to y, the derivatives of $ might 
appear through the chain rule, see also Subsection 13.21 This fact will be assumed understood, 
and it will not be mentioned anymore since it does not create any serious difficulty, or any 
difference in the estimates. 

By our construction of *$>i, all the terms multiplying powers of e less or equal than 2 reduce to 
e 2 ^fii (—A^^dhWo) + ipf'dh'Wo) = s 2 p^- ^iWq -1 dhWoipj 1 , so we are left to consider the powers 
(of e) of order 3 and higher. In the remainder of the proof, we use the symbol A 2 (e) to denote 
terms of order 1, e or e 2 : since they all generate a single term, we do not need to compute them 
separately. 

We begin by considering the terms where derivatives in y appear. Since & £ is linear in u, 
we can deal with each summand in separately. Looking at — ^det gA g (ipf (y)dhWo), second 
derivatives in y appear only in the expression — Vdet gg ab u a b, so from Lemma 1331 and Remark 
ET41 (b) we find that 

-Vdet^ A g (rf(y)d h w ) = A 2 (e) - 2e 3 (ot^) - ( n+1 H ab + Cn + iK) (tfj^dhW 

+ e 3 L 1 ip l +e A L 2 ipi, 

where L\ , L 2 are as in the statement of the lemma. 
Similarly one finds 

-Tdet^ A g (erf(y)g%(y, 0) = A 2 (e) - e^^d^ + e 3 L^ x + e 4 L 2 ^; 
-y/fetj A g (e 2 ^(y)g' 2 l (y, ()) = A 2 (e) + e 4 L 2 ^ + e 3 ^; 
-Tdet^ A 9 (e 2 w ^(y)5o (y, 0) = A 2 {e) - e 4 (d^ aVa ^)g + eV/^i + e b ^L 2 ^. 
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At this point we are left with the terms (of order e 3 and higher) which do not involve 
derivatives of ipi in y: these will appear as multiplicators of the summands in the expression of 
^i. The ones involving O^wq, g%, g% and g 3 are included in the expression e 3 Li?/>/, so it remains 
to consider e 2 mipf g^ . Recalling that ^det g = 1 + e( n H" + 0(e 2 ) (see the proof of Lemma I3T3]) . 
and expanding —pu p f £ as 



/> w% 1 + e(p-l)u% V+£ 2 (p-lK 2 w 2 + ^e 2 (p-l)(p-2)w p Q 3 wj +0(e 3 ), 
we obtain 

v / det^(l - K^V^No = M*) + e'Cn+iH^^g^l - pv?' 1 ) 

- e 3 p(p - l)ui^ 2 wi/i/VNo + eV^oV'z, 

where Lq is a multiplication operator with coefficients also satisfy (|112p . This concludes the 
proof of the lemma. ■ 

Next, using the above characterization, if 112 is a suitable linear combination of the ^/'s, we 
can estimate the scalar products of Ts e U2 (in H^ E ) with some other elements belonging to the 
subspaces H\, H2, H2 and #3, see (f88|) - (|90]) . 



Lemma 6.2 For some arbitrary real coefficients (cti)i and we consider functions u% G Hi, 
u 2 G H 2 and u 2 G H2 of the form 

00 e' s Ce~ k 

u 1 = ^2a j (p j (ey)u jt£ (\(\); « 2 = ^/3^ ; u 2 = ^ MT{ £ y)^i,e,m{Q)- 

0=0 1=0 e -ff + l 

We also let U3 G #3. Then, for 5 G (| + 7>§^ — 7) anc ^ 7 sufficiently small, we have the 
following relations 

(114) (Ts £ u 2 , m)^ = o(e 2 ) I 1 £ 



z=o 



f I II " 1. II //i 



(115) 



(116) 



(T^u 2 ,u 2 )h Se = C (l + o £ {l))—^2e 2 mPf; 



1=0 



(T^U2,U2)H Se =0(e 3 )^ (j> 2 + eV)A 2 



1=0 



Ce~ k 

a= £ -s+i 



0{£3)\\U2\\h s J\U2\\h Se ; 



(117) 



/ 1 



eV 2 + £' 



V) A 2 



z=o 
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Proof. We recall that, by Lemma EH (EH, (ED and dHTJ) there holds 



2 l + o e (l) 



£ 

i=o 



of; 



I- 1,2 1 + II 
Nlff Se = ^ ll^lWoH 



yi A 2 > 



(118) 



2 l + O e (l) 

\Mh Ss = Tk 



Ce~ k 

E 

l=£-S+l 



We show first (|114p . Since u\ is even in when we use the expression of & £ (^fi) in (| 1 13|) 
we have to consider only — 2e 3 dT^(Ei)dy ^ iji^djWQ = e 3 L 2 ipi and the errors e s L t ipi, since the 
products of all the others terms with u\ will vanish by oddness. Therefore we leave this term as 
it is, and we estimate the error terms only. So we get 

(r Ee « 2 ,«i)H s = -rrV'ayA / / %,s(|C|)<Aj(17) {^ L 2ipl + e A L^i + £ 4 /x;Li^ + e 5 mL 2 ^i) dyd(. 
j,i + 

Reasoning as in Lemma 15.71 (avoiding the scaling in e, which has been already taken care of) 
one can show that, for any integer m 



(119) 



\ |=0 / z=0 



From the Holder inequality and the last three formulas we deduce that 



4E(^ 6 (i + lwl 2 ) + e 8 lwl 3 + ^ 10 lwl 4 )A 2 



Now, from the Weyl's asymptotic formula and from the fact that 5 £ (| + 7, |/c — 7) , one finds 
that for / < e~ 5 there holds e 2 |/i/| 2 = o e (l)\ni\, that e 4 |w| 3 = o e (l) and that e 6 |^/| 4 = o e (l), so 
(fTHj) follows. 

We turn now to (|115|) . It is convenient first to evaluate some 1? norms. Writing & £ (^i) 



e 2 p^niVjQ 1 df l wotpj l + & £ (^i), and ^ = XeilCD^dhWo + ^i, from (|119p we find (I runs between 
and e- s ) 



(120) 



IE A*' IL . HE a^^-o 2 r2 < § E I 1 + - 2 + - 4 i^i 2 ) A 2 < § E A 2 ; 



L 2 



(121) 



|E^*C 2 - §E( £2 + e4 lwl 2 )A 2 * ^ 2 Ed + £ V 2 )A 2 ; 



(122) A6.(*0 2 < £ E ( £4 lw! 2 + e%t\ 2 + e 8 | W | 4 + 6 10 | W | 4 ) A 2 < ^ E rffc 

II Li o o 



50 



(123) ||5>6e(*i) [ 2 < §E( £6 |W| 2 + £ 8 |W| 4 + ^V/| 4 ) Pf < ^£(l/*| 2 + eVl 4 )tf 



Using the orthogonality of the ipi's, f)65|) and recalling the definition of C\ in Subsection (15. 2J) , 
we find 



(124) 



(r Es (tf,),*i)ir E . = £ 2 CoMj + (© 6 (*i),^^o)La + (6 E (*i), 



Multiplying by the coefficients /5's, using the Holder inequality and (|12U|) - (|123|) we get 



J> 2 + e V)A 2 £A 2 



+ E^a 2 E( 1+£ V)A 2 



Recalling the Weyl's asymptotic formula and the fact that 8 G (| + 7, |/c — 7), we obtain 
e 2 fj,f = o(/ij), eVf = o(/x;) for / < so the last formula implies (| 1 15 [) . 

To prove (I116p we notice that, by the orthogonality of the ifti's, the term of order e 2 in 
& £ (^>i), once multiplied by u 2 and integrated, vanishes identically. Therefore, from the Holder 
inequality, (fTTKj) and (fT23|) we find 



/ e -d \ 2 / 

(T S£ u 2 ,u 2 )h E£ =0( £ 3 )1 ^(/x 2 + eV)A 2 £ A 



z=o 



which is precisely (|116f) . 

It remains to prove (|117p . Using (|42p , the formulas in the proof of Lemma 13.31 and the fact 
that (linearizing ([3]) at wq) —A^(dhWo) + dhWo = pvjQ 1 dhWo, one finds 



^detg £ (-A ge ^i + ^i) = pw% $d h w + eL ^i + e 2 (L 2 ^ + uMi) + e 3 L 2 ^ 

(125) + e 4 ( W L 2 ^ + L 3 ^). 

Hence from (1113j) it follows that 

>Co /ij 



+ e 5 fiiL 2 Tpi + e 6 ^(niL 2 ipi + L 3 V0 + ©e(*i)- 

Since 1*3 is orthogonal to -ff 2 in , integrating by parts we have J s «3(— A^^+^^i/Set g £ dyd( 
for Z = 0, . . . , e" 5 . Hence from ([119]) and ([I25jl we get 



(T Ee n 2 ,«3)H Ee =0(l)||«3k Ee ( iEW + eVf + e 12 M?)/3? 



As shown before, £ /U 2 = o e (l) for I < e , so we have e \K = o(e /if), and the conclusion 
holds. ■ 



We have now the counterpart of Lemma 16.21 with u 2 replacing u 2 . 
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Lemma 6.3 For some arbitrary real coefficients and (fii)i, we consider functions u\ E H\, 
u 2 6 H 2 and u 2 £ H 2 of the form 



Ce~ 



ui = ^2aj^j(£y)u j>e (\C\); u 2 = ^2/3i^i; u 2 = ^ /3^f l (ey)^ i£)m (C). 

3=0 1=0 l=e- s +l 

Suppose also that U3 S #3. Then, for 5 G (| + 7, §fc — 7) and 7 sufficiently small, we have the 
following relations 



(126) 



(127) 



(T Ee u 2 ,u 1 ) Hs =0(e 1 -^\\u 1 



■? 1 • 



i= e -*+l 



C 



.1 Ce~ fc 



l=e-S+l 



(128) 



(T Se u 2 ,u 3 k Ee =0(e 1 ^)\\u 3 \\ H ^ 1 £ A 



i= E -<5+l 



Proof. We show first (|126p . Since u\ and «2, for any fixed y are linear combinations of 
spherical harmonics (in X) of different type, from the arguments of Subsection 14.21 it follows 
that 

/ urtCDuiMVk = 0, 



(ui,u 2 ) Hs = 0; 



so we clearly have that (Ts s ui,u 2 )h s = 0. Then (|126p follows immediately from Lemma [5TT 
To prove (|127|) , we reason as for the proof of Lemma 15.61 to find 

(129) (T Se u 2 , w) Hse = Ax + A 2 + As, 

where w G Hg e is arbitrary, and where 

Ce~ k 



Ai(w) 



-A c + (l + e 2 ^) -pu% 



m=l 



l 2 (w)=e 2 [ £ (EA((®-Wir(ey)MKI)^ 

i=£ -4 + 1 V TO= 1 I 



J=e-*+l 

As for (I95p . since — ct>j| is uniformly bounded one finds 

(130) IA2MI + |I 3 H| < Ce 2 ||u 2 || Hse ||to|| H 
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for a fixed positive constant C. Taking w = u^, by the orthogonality of the tpi's, by the fact 
that T e a Ul Vi E)m = s 2u, e v l,e,m ( see Proposition 14, 5p and by (J7U), with an integration by parts we 
have 

l=e- s +l l= £ -s + i 

From (|28p . Proposition 14.21 and Proposition 14.51 which provide estimates on <r e 3 w , e , we obtain 
(131) ii(n 2 )>^ ^ eViA 2 

for some fixed C > 0. Then (|127|) follows from (|13tjp . (j 1 3 1 j) . Lemma 15.61 and Lemma 15. II (since 
E 2 fj>i S> e 1 ' for I > and for 7 sufficiently small). 

We turn now to (|128p . By (|13U|) . taking w = u 3 , it is sufficient to estimate Ai(u 3 ) + A 3 {u 3 ). 
From r e 2 w ,uj jE)m = c e 2 w , e wj ie , m in H E 2 U E , with an integration by parts we find 



A 1 (u 3 )+A 3 (u 3 )= [ £ ^ (-A c + (l + £ 2 ^)-^ _1 ) ( E ACMMKI) fe) 



U 3 . 



From (f6"T|) and from the fact that —A^ipi = niifii + (Oi — one finds 

eWTMici) |f = -£ 2 a£cmicd y\ + e2 ((^- ®Mr^(ici) 

Therefore, integrating by parts we obtain 

(132) A 1 (u 3 ) + A 3 (u 3 ) = (U 2 ,u 3 ) Hse +A 4 {u 3 ), 

where _ 

A 4 {u 3 ) = e 2 I £ (E^, e A((«-^r(^)«z >e (KI)fe)«3, 

and where U2 = a e 2 uji ePl^^^y^l^tmiO £ #2- Now, as for U2 it is possible to prove 

that there exists a fixed C > such that 

\m\ 2 u Se <% E -Wa 2 <§ E a 2 , 

J=e-*+l l= £ - s +l 

where we used the fact that c £ 2 W; e is uniformly bounded for / < Ce~ k . Since U3 is orthogonal 
in H-£ e to H2, from Lemma 15.14 these observations and the last two formulas it follows that 

/ Ce~ k \ 5 

(U2,u 3 ) HSe = 0(s 1 ^)\\U 2 \\ H3 Ju 3 \\ H3e < Ce 1 ^ I E A 2 IMff*,. 
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1 

The arguments of the proof of Lemma [5T61 yield ^4(^3) < Ce 4 [Ylf=e- s +i 2 ll u 3ll-Hs- • Hence 
from (|129|) . (|132|) and Lemma 15. II we find that 




(Tz £ u 2 ,u 3 ) Hse = (U 2 ,u 3 ) HSe +0(e 1 7 ) I fit I ll"3||H Se 

which concludes the proof. ■ 
6.2 Applications 

In this subsection we apply the estimates in Lemmas 15. 1\ 16.21 and 16.31 to estimate the morse 
index of Ts e as e tends to zero, and to characterize the eigenfunctions of Ts e corresponding to 
resonant eigenvalues. 

From Proposition S2] we know that there exists a unique positive number a such that rfc = 0. 
If Ck is the constant given in (i27|) . we also let 

k 

(133) 9 = (JL^ 2 Vol(K). 

Then we have the following result. 

Proposition 6.4 Let O be the constant given in ()133[) . and let Ts e be the operator given in 
(|8ip . Then, as e tends to zero, the Morse index of T^ E is asymptotic to @e~ k . 



Proof. For any m € N, the m-th eigenvalue X m of Ts e , and the m-th eigenvalue X m of Tg s 
can be evaluated via the classical Rayleigh quotients 

(134) X m = inf sup — r -; A m = inf sup 



dimM m =m ueMm {u,u)h^ £ &\xaM m =m ueMm (u,u)h Se 

where M m is a vector subspace of H% £ . Choosing M m = M m to be the span of the first m 
eigenfunctions of Ts e , from the above formula for X m and from Lemma 15. II we get 

A < SUD {T ^ U)H ^ - SUD Fs e u,u)H Se +0{e 1 -">)(v,,u)H Se r Q , 

*m S SUp - - — SUp n ^ A m + <At J- 

«GM m „ e M ro (l + O^ 1 7 ))(n,n)// 5£ 

Reasoning in the same way we also find X m < X m + 0(e 1-7 ), and hence it follows that 

(135) \X m — X m \ < Ce 1-7 for all m G N and for e small, 

where C > is a fixed constant. 

Now we let N±(e) denote the number of eigenvalues X m less or equal than —e^~, and by 

~ 1 — -V 

N2{e) the number of eigenvalues X m less or equal than . From Proposition 14.61 it follows 
that N\(e) is the number of the r/; E 's which are smaller than — . Reasoning as in Corollary 
14.81 one finds that, as e tends to zero 



N 1 (e)^[ — ) Vol(K)e- 
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On the other hand, still by Proposition 14.61 we have that N 2 (e) = AT 2 + N 22 (s), where 

1—7 

A^2,i(e) is the number of rji^s which are smaller than and A^ E the number of o"; j£ 's which 

1— 7 

are smaller than . From (|27p . (|28p and Proposition 14,51 we obtain, for e small 

1 



N^-f-g) Vol(K)e- k ; 
From the last formula we deduce that also 



N 2 , 2 {e) 



k 

a n 2 



Cn- 



i.k 



Vol(K)e hi ^ Jl - k = o{e- k ) 



N 2 (e)~[^-\ Vol{K)e- k . 
Since by f)135j) the Morse index of T% e is between Aq(e) and N 2 (e), the conclusion follows. 



We can now characterize the eigenfunctions of Ts £ corresponding to eigenvalues close to zero. 

Proposition 6.5 For e sufficiently small, let A be an eigenvalue of 'Ts e such that |A| < , /or 
some if > 2, and let u £ H% e be an eigenfunction of Ts e corresponding to A with \\u\\h s = 
1. In the above notation, let u = u\ + u 2 + u%, with Ui £ Hi, i = 1,2,3. Then, if u\ = 
T l f=Q a o < t ) j{ £ y) u jfi{\C\), one has 



(136) 



r-\vjA< e ~ 



1 — 



as e — > 0. 



Proof. We show that u 2 , 113 tend to zero as e tends to zero. This clearly implies \\u— Ui\\H Se —> 
0. Once this verified, (|136p can be proved as in [30] Proposition 4.1. 

To prove that 113 tends to zero as e — > 0, we take the scalar product of the eigenvalue equation 
T-£ e u = Xu with U3. Using the above arguments (in particular Lemma l5.ip we easily find 

— +0(£ 1 " 7 )|M| //e JM| //e£ < (Tx s u,u 3 ) H v e = \(u,u 3 ) H ^ = A||u 3 ||h e • 
CC k 

This implies ||«3[|# = 0(^ 1_7 )||^||/i" Ee H^ll-H^ > an( ^ hence ||«3||jf Ee < Ce 1-7 !!^!!^^ < Ce 1_7 . 
Next we take the scalar product of the eigenvalue equation with u 2 . From Lemmas 16.21 and 
.31 we find 



frp , ^ C (l+O e (l)) \- 2 



0(1) 



Z=0 



J=0 / \ l=e- s +l 



+ 



c 



_1 Ce" 



E £2 wA 2 - 



Since e 5 /^ 2 + e 7 [i\ = o e (l)|^| for / < e~ s and e = o(e 2 /i/) for / > e~ s (recall that 5 G 
(| + 7, k — 7)), it follows that 



(137) 



(T^u 2 ,u 2 ) Hs£ > C- 1 -^ E £2 WA 2 



=0 
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for a fixed positive constant C. Finally, still from Lemmas I6.2H6.31 from the fact that e 4 |/x;| + 
£ 6 |w| 3 = °e(l) f° r ^ — £ ~ S an d £ 2 ~ 27 = o(e 2 fn) S> 1 for I > (taking 7 sufficiently small) we 
have also that 



(138) (T Se n 2 ,ui + n 3 )// Ee = o e (l)(||wi||H Be + ||u 3 ||ff Be ) ( ^ ^N^J 



1=0 



From (11371) and (j!38j) and the fact that Ts £ is self-adjoint we deduce that 

c- 1 



_ 1 

Ce- k I Ce~ k \ 2 



^2 e 2 fj, t /3f + o e (l) [ e 2 |^|A 2 ] (ll«illH Se + ||«3||ftj«) < {Tt, £ u,u 2 )h^ = X(u,u 2 )h Ss 

< Ce^WuWh^ \\u 2 \\h Ss - 



(=0 \ 1=0 



Also, from Lemma \5A\ testing the eigenvalue equation on Yli<i Pi^h where Iq is the biggest 
integer such that pn Q < 0, one finds 



e 

l<l 



The last two formulas imply that ^ A 2 = °e(-Q> namely that [|ii2 [|^r Se t en ds to zero as e 

tends to zero. This concludes the proof. ■ 



6.3 Proof of Theorem fTTTI 

Once Propositions 16.41 and 16.51 have been established, the proof goes as in [39], Section 8 (see 
also [38] Section 5) and therefore we will limit ourselves to sketch the main steps. 

First of all, using Kato's theorem, see [30], pag. 445, one can prove that the eigenvalues of 
are differentiable with respect to e, and if A is such an eigenvalue, then there holds 

(139) — = {eigenvalues of Q\} , 
where Q\ : H\ x H\ — > R is the quadratic form given by 

(140) Q x (u,v) = (1- A)| Vu-Vv-p(p-l)J^ uvu% 2 (^^{e-). 

Here H\ C H% £ stands for the eigenspace of Ts e corresponding to A and the function uj t£ : Q — > R 
is defined by the scaling ui )£ (x) = n/ i£ (ex), where ui jE is as in Section [3j Notice that, since A 
might have multiplicity bigger than 1, when we vary e this eigenvalue can split into a multiplet, 
which is allowed by formula (|139H . 

Taking A as in Proposition 16.51 we can apply (|139p . and evaluate the quadratic form in (|140p 
on the couples of eigenfunctions in H\, which are characterized by (|136p . Reasoning as in [38], 
Proposition 5.1 one can prove the following result. 
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Proposition 6.6 Let A be as in Proposition 1 6. 51 Then for e small one has 

<9A 1 .— ... 

where F is a positive constant depending on N, k and p. 

Now we are in position to prove the following proposition, which states the invertibility of T% e 
for suitable values of e. 

Proposition 6.7 For a suitable sequence Sj — ► 0, the operator J £ (ui £ ) : -fr 1 (f2 e ) — ► is 



invertible and the inverse operator satisfies 



< . Fh <n > f or all 3 £ 

m(n s .) ~ nun{ £ f, £ J}' 



Proof. From Proposition 16.41 we have that, letting iV e denote the Morse index of Ts e , there 
/-\- 

holds N £ ~ 2 Vol(K)e- k . For I G N, let £/ = 2~ l . Then we have 



k . fc 



(141) iV £;+i _ at £; ~ ^ 2 y /(K)(2^ +1 ) _ 2 M ) ~ ^ 2 Fo/(^)(2 fc - l) e 7*. 

By Proposition 16.61 the eigenvalues A of Ts e with |A| < £ ? are strictly monotone functions of e 
so by the last equation the number of eigenvalues which cross 0, when e decreases from e\ to 
£/+i, is of order e7 k ■ Now we define 

Ai = {ee (e I+ i,£j) : kerT S£ /0}; Bj = (e J+ i,ej) \ A- 

By Proposition 16.61 and (|14ip we deduce that card(j4/) < Ce^ k , and hence there exists an interval 
(ai,bi) such that 

(142) k«S* ft-^C-^l^-.^. 

From Proposition 16.61 then it follows that every eigenvalue of Ts a +b in absolute value is bigger 

— 2 — 

than C _1 min{e fc , £ f } for some C > 0. By Lemma 15.21 then the same is true for the eigenvalues 
of J"(ui >£ ) so the conclusion follows taking ej = aj ^" 6j . ■ 



Remark 6.8 The arguments in the proof of Proposition \6.5\ can be easily adapted to the case 
in which |A| < C~ 1 e 2 with C is sufficiently large. Therefore the result of Proposition pT7| can be 



improved to 



< ■ Fk in , for all j G N. 



Below, || • || denotes the standard norm of H (fl £ ). For the values of e such that J £ (uj s ) is 
invertible, it is suf ficient to apply the contraction mapping theorem. Writing £ = £,-, we find a 



solution u £ of (P^) in the form u £ = Uj £ + w, with w G H l (Q, £ ) small in norm. Since J £ (uj j£ ) is 



invertible we have that J' £ (u) = if and only if w = — (J £ (uj j£ )) 1 [J' £ {ui )£ ) + G(w)], where 

G(w) = J' £ (u I;£ +w)- J' e {u I>£ ) - J'J(u I>£ )[w}. 
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Note that 



G(w)[v] 



(ui jE + w) p -u p j £ - pu p j 



p-1 



w 



v. 



v£H\n £ ). 



Reasoning as in the last section of |40| . we find the following estimates, which are based on 
elementary inequalities 



(143) 

(144) 

\\G(W) 



\\GHW < 



C\\w\\ p for p < 2, 
C||w|| 2 for p > 2; 



Ml < 1; 



G(w 2 )\\ < 



\ w i\ 



\w 2 \\ < 1. 



C (IHf- 1 + IKf" 1 ) IK - w 2 \\ p<2, 
C (||u>i|| + \\w2W) \\wi — w 2 \\ p > 2; 

Defining F £ : H 1 ^) -» H 1 ^) as 

F e (w) = - {Ji'^))' 1 [J' s (ui >£ ) + G(w)] , w e H 1 ^ 

we will show that F E is a contraction in some closed ball of i^ 1 (0 £ ). From (|40p . Proposition [6T 
(with Remark ESI) and p33" ]) - ([H4]l we get 



(145) 



\FJw)\\ < I 




1+1- 



1+1- 



2 + \\w\ 



2 + \\w\ 



for p <2, 
for p > 2; 



(146) 



Ce^+V (llu-illf- 1 + ||w 2 
v Ce-( fc+1 )(|H|| + ||«; 2 || 

Now we choose integers d and k such that 

for p < 2, 



IP- 1 ! 



\W\ — W2\ 



\W\ — W2\ 



P<2, 
p>2; 



w\\ < 1; 



\Wi\ 



\uu2W < 1. 



(147) 
and we set 




for p > 2; 



I > d 



Iwll < £' 



'}■ 



From (|145p - (|146p we deduce that is a contraction in S for e small, so the existence of a critical 
point u £ of J e near follows. All the properties listed in Theorem ll.il including the positivity 
of the solutions, follow from the construction of u/ i£ and standard arguments. As in [3D], when p 
is supercritical one can use truncations and L°° estimates to apply the above argument working 
in the function space H l (Q £ ) n L°°(Q £ ). 

Remark 6.9 With the arguments given in Section we could obtain sharp estimates on the 
Morse index of and on the eigenf unctions corresponding to resonant eigenvalues. In par- 
ticular about the latter we showed that the components in H2,H^ are small, and that in Hi the 
Fourier modes are localized near some precise frequencies. This allowed us to prove Proposition 



6. 7| using Kato 's theorem. 

Even if we did not work the computations out, it seems it should be possible to give a more 
rough characterization of these eigenf unctions (in particular on the H 2 component) and to prove 
a (non sharp) estimate on the derivatives of the eigenvalues, still obtaining invertibility. This 
might slightly simplify the proof of existence, although most of the delicate estimates will be 
shifted from the analysis of Ts e to that of the quadratic form Q\ defined in (|140() . 
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